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ABSTRACT 

Aims. We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the 
influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. 
Methods. Using the IRAM 30 m telescope, we detected HCN, HCO+, '-CO, and ''CO in six GMCs along the major axis of M 33 at 
a resolution of ~ 1 14 pc and out to a radial distance of 3.4 kpc. Optical, far-infrared, and submillimeter data from Herschel and other 
observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of 
photon-dominated regions, one for solar and one for M 33-type subsolar metallicity. 

Results. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large 
Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The '"CO/'^CO line ratio varies between 
9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed 
HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios 
observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, 
while HCO^ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO^/HCN 
and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of 
moderate densities of n 3 10^-3 lO'' cm"', while the FUV field strength only has a small effect on the modeled line ratios. The line 
ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor 
role in influencing these line ratios. 

Key words. Galaxies: individual: M33 - Galaxies: ISM - ISM: molecules - ISM: clouds - ISM: photon-dominated region (PDR) 



1. Introduction 

Owing to their large dipole moments, even the rotational ground 
state transitions of HCN and HCO^ trace dense molecular gas 



Send offprint requests to: C. Buchbender, e-mail: buchbendOiram. es 
* Based on observations with the IRAM 30m telescope, Herschel, 
and other observatories. IRAM is supported by CNRS/INSU (France), 
the MPG (Germany), and the IGN (Spain). Herschel is an ESA 
space observatory with science instruments provided by European-led 
Principal Investigator consortia and with important participation from 
NASA. 

** FITS files of the presented spectra of the ground-state tran- 
sitions of HCN, HCO+, '-CO and "CO are available at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsw eb.u-strasbg.fr /cgi-bin/qcat?J/A-HA/| 



with densities in excess of ~ 10 cm . Because stars condense 
out of dense cores of giant molecular clouds (GMCs), both 
molecules are promising tracers of star formation (SF) and the 
star fo rm ation rate ( 'SFR) . A series of paper s (IGao & Solomor 



2004allbl: IWuetal 



2005 



Gao et al 1 120071 



ijGao & Solomod 
iBaan et al.. ,20081 



Gracia-Carpio et alj|2008t IWu et al.ll2010t iGarcfa-Burillo et all 
20121: lLiu&Gadl2012h have recently investigated the corre- 



lation of HCN (and partly HCO+) with far-infrared (FIR) lu- 
minosities (Lfir) in galactic GMCs, centers of nearby galax- 
ies, and (ultra-)luminous galaxies (LIRGs/ULIRGs), showing 
that HCN is indeed a good tracer of SF and tightly correlated 
with LpiR. Th ere are, however, findin gs that complicate this pic- 
ture (see e.g. lCostaghola et"al 1 1201 Ih . In contrast to CO, which 
traces the bulk of the molecular gas-phase carbon, HCN and 
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Fig. 1: SPIRE 250 |^m map of M 33 dXilouris et al.ll20T2l) . The 
rectangle delineates the 2' X 40' wide strip along the major axis 
shown in Fig. 2. Crosses mark the positions of the observed 
GMCs. 



scope. M 33 is a spiral galaxy with Hubble type SA(s)cd lo- 
cated at a distance of only 840 kpc (Table [T] and Fig. [T]i. It is 
the third largest member of the Local Group (after M31 and 
the Milky Way). Observations of small-scale structures in M 33 
do not suffer from distance ambiguities as galactic observations 
do. Its small distance allows us to obtain a spatial resolution of 
~ 114pc (i.e. 28") at a frequency of 89 GHz (i.e. 3.3 mm) with 
the 30 m telescope. M33 is seen at an intermediate inclination 
of ; - 56°, yielding a short line-of-sight depth, which allows 
us to study individual cloud complexes. It is roughly ten times 
less massive than the Milky Way, and its overall me tallicity is 
12-t-lo g 0/H=8.27, subsolar by about a factor two jMagrini et al.l 
2010h . Therefore M33 is particularly interesting to compare 
with the Milky Way, but also with th e Large Magellanic Cloud 
that has a metallicity similar to M 33 ( iHunter et alj|2007 ). 

Using the IRAM 30 m telescope. iRosolowskv et al I ( 1201 Ih 
(hereafter RPGll) observed four massive GMCs of more than 
3 1O^M0 in M33, searching for the ground-state transition 
HCN. They detected HCN in only two of the GMCs. The ob- 
served GMCs are under-luminous in HCN by factors between 
two to seven relative to their CO emission when compared to 
averaged values in the Milky Way. 

Here, we present new, deep observations of the ground-state 
transitions of HCN, HCO+, '^CO, and CO towards seven GMCs 
in M33 including three of the clouds observed by RPGIL All 
four tracers are detected in six of the GMCs. The relative weak- 
ness of HCN emission is confirmed and interpreted using models 
of photon-dominated regions (PDRs). To better characterize the 
observed GMCs, we estimated their star formation rate, total in- 
frared luminosities, and FUV fields using a large ancillary data 
set compiled in th e framework of the Herschel open time key 
project HerMBBes dKramer et al.ll2010l) . Figure|2]shows a subset 
of this data set. 



HCO^ are minor species. Their abundances are therefore more 
strongly influenced by the d etails of the chemical network (e.g. 
iLopez-Sepulcre et alj|2010b . HCO^ is linked via ion-molecule 
reactions to the ionization equilibrium. Its collisional cross- 
section is close to a factor 10 larger than that of HCN, which is 
linked to the hydrocarbon chemistry and the amount of nitrogen 
in the gas phase. Elemental depletion in low-metallicity environ- 
ments may therefore have a strong effect on its abundance. 

Most extragalactic observations of HCN and HCO^ have so 
far been restricted to the nuclei of galaxies or their integrated 
fluxes. Exceptions are a study of HCN and HCO^ emission 
in the disk of M3I by iBrouillet et alj (l2005h. LMC o bserva- 
tions ( IChin et al.ll 19961 1 19971 1 1 998'; 'Heikkila et al.' 1999), HCN 
maps of seven Seyfert galaxies by Curran et al. (2001), and also 
HCN mapping along the major axis in 12 nearby galaxies by 
iGao & Solomoiil (l2004albl) (hereafter GS04a,b). We also recom- 
mend the s t udies of H CN and CO and their ratios in M5 1 by 
iKuno et alj ( Il995h and lSchinnerer et al.l ( |2010|) . The intersteflar 
medium (ISM) of nuclear regions are often subject to particu- 
larly strong heating sources because they are often dominated 
by starbursts with intense UV fields heating the gas or active 
galactic nuclei (AGNs) with strong X-ray emission. Indeed, the 
HCN to HCO^ line intensity ratios are found to be systemati- 
cally higher in AGN-dominated regions, such as in the central 
part of NG C 1068, and low in pure starburst environments, as 
in M82 (e^. lKohno et al.ll2003l:llmanishi et alJl2009l:lKrips et aP 
[2008, 20ll. 

We have targeted seven GMCs along the major axis of M 33 
out to a radial distance of 4.6 kpc, using the IRAM 30 m tele- 



Table 1 : Basic properties of M 33 





M33 


References 


RA(2000) 


01:33:51.02 




DEC(2000) 


30:39:36.7 




Type 


SA(s)cd 


1 


Distance [kpc] 


840 


2 


11" (30m @ 230 GHz) equal to 


45 pc 




21" (30m @ 115 GHz) equal to 


86 pc 




28" (30m @ 89 GHz) equal to 


114pc 




LSR velocity [km s '] 


-180 




Position Angle [deg] 


22.5 


3 


Inclination [deg] 


56 


4 


Ris 


30.8' or 7.5 kpc 





Referenc es. (1) Ide Vaucouleurs et alj jl991h: (2) iGalleti et al.l i200§\ 
Freedma n^etal] 71991 ): (3) P aturel et all (l2003h : (4) .Reean & Vogel 
a994l) : IZaritskv et alj ( 119890 . 



2. IRAM 30 m observations 

We used the IRAM 30 m telescope to perform single-pointed ob- 
servations of the ground-state transitions of HCN, HCO^, '~CO, 
and '^CO towards seven GMCs in M 33. Observations were car- 
ried out between December 2008 and July 2012, comprising a 
total of 109 hours of observing time. In 2008, we used the now 
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Fig. 2: Observed positions towards seven giant molecular clouds (GMCs) within a 2' x 40' strip along the major axis of M33. 
The strip extends from 10' south of the galactic center to 33.3' north. The center of the strip is at 01:34:11.8 +30:50:23.4 (J2000). 
Circles indicate the 30 m b eam size of 28" at 90 GHz. Panels show fro m top to bottom: integr ated intensities of Ha emission 
dHoopes & WalterbosI I2OO0I) 24; 70 ixm emission observed with Spitzer jTabatabaei et al.l |2007|): continuum emission between 
100 M - m and 500 |j.m obser ved with PACS andSPIRE in the framework of the HerMB 3es pro gram dKramer e t al. B oquien et alj 
I2OI lUXilouris et al.ll2012h : '^CO 2-1 30 m observation and Hi VLA data, both taken from lOratier et al.1 (I201Q) . All data are shown 
at their original resolutions. 



decommissioned A 100 and BlOO receivers that had a bandwidth 
of 500 MHz and the IMHz filterbank, to observe each of the four 
Unes individually. 

The bulk of the observations were carried out in 2009 em- 
ploying the new eight-mixer receiver EMIR and its instanta- 
neous bandwidth of 16 GHz in each polarization, connected 
to the wide-band WILMA autocorrelator backend with 2 MHz 
spectral resolution. This setup allowed simultaneous observation 
of HCN with HCO+ and '^CO with '^CO. One advantage of the 
simultaneous observations is that the relative intensity calibra- 
tion of the lines is very accurate. The observations were carried 
out in wobbler switching mode using the maximum available 
throw of + 120 "and a switching frequency of two seconds. This 
mode ensures more stable baselines than the position-switching 
mode. However, the velocity resolution of about 6kms ' in the 
3 mm band only barely resolves t he spectral lines of M 33, which 
are typically 10-15 km s ' wide (iGratier et al.ll20ldl) . The beam 
sizes are 21" at 115 GHz and 28" at 89.5 GHz, corresponding 
to a spatial resolution of 86 pc and 1 14 pc, respectively, in M 33 
(cf. Table[B. 

The observations of '"CO and '^CO were repeated in June 
and July 2012 using position-switching and the FTS spectrome- 
ters with an off position outside of the disk of M33 to exclude 
the possibility of self-chopping effects in the spectra. The latter 
were present in some of the earlier wobbler- switched '"CO spec- 
tra. Due to the high critical density of the dense gas tracers HCN 
and HCO^, as well as the observed velocity gradient along the 
major axis of M33 (see Table |2|, we reckon that self-chopping 



is less probable for these lines. Please note that the '^CO data of 
GMCl, GMC91, and GMC26 are taken from RPGll who also 
used position-switching. 

All data were reduced using the GILDAS Q software pack- 
age. Each scan was inspected and scans with poor baselines 
or unreasonably high rms values were rejected. Before averag- 
ing, linear baselines were fitted and removed. The data were 
regridded to a common velocity resolution. Spectra were con- 
verted from the to the r,nb scale by multiplying with the ra- 
tio of forward efficiency (Fgif = 95%) to main beam efficiency 
(Beif - 81%), taken to be constant for the observed 3 mm lines. 
The reduced spectra are shown in Fig.|3] 

Integrated intensities were extracted from the spectra on 
a scale by summing all channels inside a velocity range 
around each particular line. The velocity range was determined 
by eye for each position from the full width to zero intensity 
(FWZI) of the '^CO 1-0 Une and is marked in Fig. [3] We deter- 
mined cr uncertainties of the integrated intensities by measuring 
the baseline rms (7"™) in a 300 km s ' window centered on the 
particular line using the corresponding '^CO 1-0 FWZI as base- 
line window and using cr - T^^^ y/N Av^^s with the number of 
channels and the channel width Av,es- In case the 1 cr value is 
higher than the typical 10 % calibration error of the IRAM 30 m, 
the former is used to estimate the observational error for the fol- 
lowing analysis. If the integrated intensities are lower than 3 cr, 
we use this value as an upper limit. Table|2]lists the observed in- 



' http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 3: Spectra of the ground-state transitions of HCN, HCO^, '^CO, and '^CO at the positions of seven GMCs along the major 
axis of M33 (cf. Fig. 1, 2). '^CO and '^CO have been observed in position switching; HCN and HCO^ with wobbler switching. 
The spectra are shown on a main beam brightness temperature scale (Tnib)- The velocity resolution is given by the spectrometer 
with the lowest resolution, i.e. WILMA, and is 5.4 km s ' in case of ''CO and '^CO, and 6.7 km s ' for HCN and HCO^. Center 
velocities are listed in Table |2] The local standard of rest (Isr) velocity displayed covers 300 km s"'. The same velocity range is 
used to determine the baseline (red lines), excluding the line windows determined from '^CO 1-0 (dotted lines), cf. Fig. lA.ll 



tensities, intensity ratios, and further ancillary data. For details 
on the latter see AppendixICl Error estimates are given in paren- 
theses after the integrated intensities in Table |2] 



3. GMCs: Selection of positions and properties 

Motivated by the HerM33es project, the GMCs were selected 
to lie within a 2' x 40' wide strip along the major axis of M 33 
shown in Figs.[T]and|2]at a range of galacto-centric distances of 
up to 4.6 kpc. Three of the GMCs (GMCl, GMC26, GMC91) 
belong to the sample of CO-bright clouds studied by RPGl 1 in 
search of HCN emission. We added four other GMCs (no6, no3, 
nol, no2) to increase the range of studied galacto-centric radii, 
as well as physical conditions. 

Table |2] lists their observed properties and Appendix |C] de- 
scribes in detail how they were derived. The masses of the 
molecular gas traced by CO, calculated using Xco-factors de- 
rived individually for every cloud as a function of integrated 
CO 1-0 intensities and total IR luminosity (cf. Appendix IC.4b . 
vary by a factor 130 between 0.1 10^ (GMC no2) and 13 10^ Mq 
(GMCl). The SFRs vary by more than a factor 50 and the 
far ultraviolet (FUV) field strengths by a factor larger than 20. 
The GMC near the nucleus, GMCl, is the most massive and 
shows the strongest SFR, as well as the highest FUV flux, while 
GMC no2 at 4.6 kpc radial distance is the least massive in molec- 
ular mass and shows only weak activity. 

Individual areas of the strip have been mapped in [C ii] and 
other FIR lines in the framework of the HerM3 3es project, which 



will yield additional insight into the properties of the ISM of 
M 33. In the first papers on fC ii], we focused on the Hii region 
BCLMP302 (iMookeriea et al.11201 ll) which is associated with 
GMCno3, and on BCLMP691 (Braine et al., in subm.), which 
lie s near GMC nol. 

iGratier et all (1201 2|) identified over three hundred '^CO 2-1 
clumps in M 33 and present a detailed view of each individual 
clump in Ha, 8 i.im, 24 [a.m, and FUV, together with the corre- 
sponding HI and '^CO 2-1 spectra and further complementary 
data. The seven GMCs discussed here are among the identified 
clumps. In Table|2]we give the corresponding clump numbers. 



4. Observed line ratios 

4.1. Spectra at individual positions: HCO+ and HCN 

HCO^ is detected at six positions with 6 to 12 mK peak tem- 
peratures and with good signal-to-noise ratios of at least seven; 
position no2 has not been detected. HCN emission is detected 
at the same six positions with signal-to-noise ratios of four and 
better; position no2 is but tentatively detected at a signal-to-noise 
ratio of 3.5. The HCO^/HCN ratio of line integrated intensities 
of positions where both molecules are detected varies between 
1.1 and 2.5 (Table|2]i. Below, we compare the observed ratios in 
detail with ratios found in the Milky Way and in other galaxies. 
Although the integrated intensities we find for GMC26, GMCl, 
and GMC91 differ up to a factor of two from the values and up- 
per limits given in RPGl 1 for the same positions, they are con- 
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Table 2: Observed intensities and complementary data 
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"4-. "4- (^U. 1 ) 


C A /A OA 

J. 7 VU.Z) 


1.4 1,U.UJ 


9 7 /-A 1 A 
Z. / l^U. 1 J 


l.J vU.UJ 


V.I (U.i) 


A a /A AA 

U.J (U.U) 




3.5 (0.7) 


2.4 (0.3) 


1.6 (0.3) 


3.0(0.5) 


1.3 (0.3) 


2.0 (0.4) 


0.8 (0.2) 


^tir/^hco+ [10 ] 


1.4 (0.1) 


2.2 (0.2) 


1.4 (0.1) 


1.5 (0.2) 


0.9 (0.1) 


1.5 (0.2) 


> 1.6 


SFR [Mq Gyr-' pc'-] 


35.9 (4.3) 


65.0 (7.8) 


6.6 (0.8) 


13.7(1.8) 


4.0 (0.6) 


12.2(1.7) 


1.2 (0.1) 




5.1 


6.9 


1.6 


3.2 


1.5 


2.0 


0.3 


Mhi [10^^ Mq] 


9.4 (1.4) 


5.8 (0.9) 


4.1 (0.6) 


8.0(1.2) 


8.8(1.3) 


5.2 (0.8) 


4.8 (0.7) 


Mh, [10= Mq] 


9.6(1.0) 


12.7 (1.3) 


3.0 (0.3) 


7.9 (0.8) 


8.2 (0.8) 


2.1 (0.2) 


0.1 (0.0) 


Go 


37.3 (1.2) 


50.7 (1.6) 


11.6 (0.4) 


23.5 (0.8) 


11.3 (0.4) 


14.4 (0.5) 


2.5 (0.1) 


Av 


6.3 (0.6) 


6.1 (0.5) 


2.3 (0.2) 


5.3 (0.5) 


5.7 (0.5) 


2.4 (0.3) 


1.6 (0.2) 



Notes. Top panel: line intensities are on the T^jb-scale and on their original resolutions: 12" for '^CO 2-1, 24" for '*C0 and "CO 1-0 and 28" 
for HCN and HCO^ 1-0. Bottom Panel: line ratios and complementary data are o n a common resolution of 28". See Appendix |C] for details. 

'^CO 2-1 clump numbers fr om Gratier et al. (2013); "'' iRosolowskv et al.l l [20TIh ; <'> FWHMs of Gaussian fits to the high resolution '^CO 1-0 
spectra (Fig. IA.lt : ^'^' iGratieret al.. 1201 ft) : Baseline rms of HCO* spectra at a velocity resolution of 6.7 km s"'. 



sistent within 3cr of the baseline rms of the observations from 
RPGll. We attribute the discrepancies to baseline problems of 
the RPGll data. 



4.2. Spectra at individual positions: ^^CO and '■'CO 

Emission from '-CO and '^CO is detected at all seven positions, 
though varying by a factor of more than 15 between GMCno2 
and GMC91. The ratio of line integrated ''CO vs. '''CO inten- 
sities varies between 9 for GMCnol and 15 for GMCno6. In 
Fig. lA. II we show the '^CO and '^CO spectra at a high resolu- 
tion of 1 kms '. The resolved line shapes are Gaussian and the 
corresponding FWHMs are given in Table 2. 



4.3. Stacked spectra 

Stacking of all spectra allows improvement on the signal-to- 
noise ratios and detection of faint lines. Figure|4a] shows the 
stacked spectrum of all data taken near 89 GHz. It was created 
by shifting all spectra in frequency such that the emission lines 
align in frequency. Individual spectra were weighted by integra- 
tion time. In addition to the lines of HCN and HCO^, the re- 
sulting spectrum shows detections of CCH and HNC 1-0. The 
average baseline rms is 0.27 mK at 5.4km s"' resolution. HOC^ 
is tentatively detected with an upper limit of 19 mKkms"', re- 
sulting in a lower Hmit to the HCO+/HOC+ ratio of 5.8. 



The stacked spectrum centered on 112 GHz (Fig.l4b]i does 
not show additional detections other than '^CO and '^CO even 
after additional smoothing of the velocity resolution. Table |3] 
lists the integrated intensities and upper limits of all detected 
transitions in the stacked spectrum, as well as their correspond- 
ing LTE column densities and abundances, derived as explained 
in Appendix Id] 

4.4. Comparison witli other sources 
4.4.1. HCO+ZCOvs. HCN/CO 

For the comparison of different tracers, all data were convolved 
to the same resolution of 28". We account for the different 
intrinsic beam sizes of the CO, HCN, and HCO^ 1-0 obser- 
vations by multiplying with beam filling factors determined 
from the '-CO 2-1 map dGratier et al.ll201()l) . Figure|5]compares 
the HCO+/CO vs. HCN/CO line intensities ratios observed in 
M 33, with those observ ed at nine positions i n the disk of the 
Andromeda galaxy M 3 1 (iBrouillet et al.ll2005L hereafter BR05). 
M 31 lies at a similar distance as M 33 of 780 kpc and had been 
observed with the 30 m telescope as well. Therefore, both studies 
obtain about the same spatial resolution of ~ 1 14 pc. 

For M33, we find HCN/CO ratios in the range of 
0.4% -2.9% (mean: 1.5 + 0.8%) and HCO+/CO ratios in 
0.6% -3.5% (mean: 1.9 + 1.0%). BR05 finds comparable 
values in the spiral arms of M31: HCN/CO 0.75% -2.8% 
(mean: 1.7 + 0.5%) and HCOVCO 1.1% -3.9% (mean: 
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Table 3: LTE column densities from the stacked spectra. 





C2H 


HCN 


HCO+ 


HOC+ 


HNC 


''CO(l-O) 


'^CO(l-O) 




/ [rnKkms"'] 


25.1 


81.0 


110.0 


< 18.8 


13.7 


522.0 


6280.0 




A'[x]* 


5.10e+12 


8.03e+ll 


2.11e-Hll 


< 7.10e-hl0 


4.18e+10 


1.80e+15 


8.55e-Hl5 


8.47e+20 


A'[x]/A'h/' 


-8.40 


-9.20 


-9.78 


< -10.25 


-10.48 


-5.85 


-5.17 


1.00 



Notes. deduced from '^CO, see Appendix |DJ *** column density; logarithmic values. 




8BOOO B9000 90000 

Rest Frequency (MHz) 



(a) 




112000 
Rest Frequency (MHz) 



(b) 

Fig. 4: (a) Stacked spectrum of all data taken in the frequency 
range between 87.2 and 90.8 GHz. The red dashed line denotes 
the 3 cr value average over the entire baseline. The HCN and 
HCO^ lines are not shown up to their maximum peak temper- 
ature, (b) Stacked spectrum of the wobbler switched data taken 
in the frequency range between 108.1 and 115.5 GHz. The av- 
erage 3 0" value is shown as red dashed line. The baseline noise 
increases with frequency because of the increasing atmospheric 
opacity. C'^O, C'^O, and CN are marked but not detected. The 
'"CO and '^CO lines are not shown up to their maximum peak 
temperature. 



2.0+0.7%). A linear least squares fit to the M33 data results in 
HCO+/CO = (1.14±0.15%)HCN/CO + (0.18 + 0.14%). This 
is consistent within errors to the fit results obtained by BR05 for 
the M 31 data: HCOVCO = 1.07%HCN/CO + 0.23%. We 
excluded positions with upper limits from the fit, i.e. position 
GMCno2. 

In the Milky Way in the solar neighborhood values of 
HCN/CO are found between 0.7% - 1.9% (mean: 1.4 ± 
2%), while the Gala ctic plane hosts on average 2.6 ± 0.8 % 
dHelfer & Blitzlll997l) . HCO+/C O values in the Galac tic center 
range between 0.9 % and 7.6 % jRiquelme et al.ll2010l). 

The HCN/CO ratios found in the LMC bv lChin et all (Il998h 
and iHei kkila et al. ( 1999) lie between 3 % and 6 %, and are thus 
higher than any value found in M33, M31, and also M51 
where HCN/CO - l.l%-2% in the spiral arms are reported 



o 
u 



o 

U 2 
X 



— /(HCN)=/(HCO+) 

— FitM33 

— ■ FitM31 (BROS) 

□ Galaxies GA04 

□ MWDiskHE97 
* M31BR05 

M33 thiswork 



1 2 3 4 5 

/(HCN)//(CO) [%] 

Fig. 5: Ratios of integrated intensities HCO+/CO vs. HCN/CO 
for M33 (red points: this work) and for M31 (green points: 
BR05). Upper/lower limits are denoted by arrows. Linear least 
squares fits to data from M 33 (red solid line) and M 31 (BR05, 
black dashed) are shown. Both fits exclude points with upper 
limits. The solid black line shows the angle bisector where 
/(HCO^) - /(HCN). The gray shaded areas display the range 
of the HCN/Cq foun d in the disk of the Milky Way (MW) by 
(iHelfer&Blitjl 19971 HE97) (darker gray) and in a sample of 
normal spiral galaxies (GA04b) (lighter gray). 



(iKuno et al.lll995h . Unlike GS04b in a sample of normal galax- 
ies, we do not find a systematic change in the HCN/CO ratio 
between regions in the center of M 33, i.e. the inner ~ 1 kpc 
(here GMC 1 , GMC26) and regions at greater galacto-centric dis- 
tances (cf. Tabled. GS04b reports that HCN/CO drops from ~ 
10% in the centers of normal galaxies to ~ 1.5% - 3% in their 
disks ^ 4 kpc. In ULIRGs and AGNs the ratios may reach global 
HCN/CO values as high as 25% (GS04b and references therein). 
GS04b attribute these high ratios to the presence of starbursts 
and argue that HCN/CO may serve as a starburst indicator. 

4.4.2. HCO+/HCN vs. HNC/HCO+ 

The HCO^/HCN ratios observed in M33 vary between 1.1 
and 2.5, while the upper limits derived for the HNC/HCO^ ra- 
tios vary between 0.2 and 0.5 (Fig. |6]l. The upper limit of the 
HNC/HCN ratio is at maximum 0.7 (GMC91), while the stacked 
spectrum where HNC has been detected shows a HNC/HCN ra- 
tio of 0.17 (Fig.©. 
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z 
u 
a: 



o 
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X 



0.4 



0.2 



0.0 



- -0.2 



-0.4 






^ Arp220 


♦ LMC (CH97, CH98) 

□ LIRGSAJLIRGS (BAIO) 
X LIRGSAJLIRGS (CO 11 ) 

□ Range M3 1 (BROS) 
M33GMCS 

4h M33 stacked 


B-^ □ f. 

□ \" ^\ 

%r % ^% 
\ % % 


-1.0 -0.5 0.0 0.5 



logHNC(l-0)/HCO+(l-0) 

Fig. 6: Comparison of integrated intensities HCO^/HCN vs. 
HNC/HCO^ in M33 (red filled circles and square; arrows 
indicate upper/lower lim its) with values found in the LMC 
(IChin etalJl 19971 Il998h (CH97, CH98; bl ue diamonds) and 
in luminous infrared galaxies compiled by Baan et al.' (l2008h 
(BA08; open symbols) and by [Costagliola et al. (2011) (COll; 
crosses). The dotted, dashed, and dot-dashed lines shows 
HNC/HCN = 1, 0.4, and 0.17 (stacked value of M33), respec- 
tively. The gray shaded area shows the range of the observed 
HCOVHCN ratios in M 31 by BROS. 



These ratios are c ompared wi t h the ratio s found in lumi- 
nous i nfrared nuclei bv lBaan et alJ (|2008|) and Costagliola et alJ 
(1201 Ih and i n Hn regions of the LMC bv lChin et alJ(ll997lll998h 
(cf. Fig.3a in lBaan et all (l2010h ). The HCO+/HCN ratios are also 
compared to the range found in the disk of M 3 1 by BROS (cf. 
Fig.©. 

The HCO+/HCN ratios, found in M 33 in the six GMCs with 
clear detections, lie at the upper end of the distribution of val- 
ues found in LIRGs. While the starburst galaxy M 82 exhibits a 
higher ratio than the AGN NGC 1068 (1.6 vs. 0.9), these ratios 
lie within the scatter of values and their errors observed in the 
disk ofM33. 

The HCO^/HCN ratios in M 33 lie in the overlap region be- 
tween the ones found in M 31 and those found in the LMC. We 
find neither ratios as high as in the LMC, i.e. 3.S, nor ratios as 
low as in M 31, i.e. O.Sl. Interestingly, all detected regions in the 
LMC are situated in the same parameter space as those detected 
inM33. 

The detection of HNC in the stacked spectrum allows us 
to derive an average HCO+/HNC ratio of 7.8 (Table[3l) for the 
GMCs observed in M 33. This ratio lies at the very high end of 
the range of values found in any of the other samples plotted in 
Fig.lSl 

More remarkably, the HCN/HNC value of S.9 from the 
stacked spectra is higher than any from the other samples we 
compare with in Fig.|6] Furthermore, it is higher than ratios 
observed over the surf ace of IC 342 which are only ~ 1 - 2 
eier & Turnej llOOSh . higher than the ratios in a range of 



galaxies found bv lHuettemeister et al.l (ll99Sh of < 4, higher than 
the typical ratios of 1 observed in starburst and Seyfert galaxies 
(e.g. Aalto et al. 2002), and also higher than ratios of 1-3 o b- 
served in Galactic molecular complexes dWootten et al.lll978h . 

This extraordinary high ratio indicates that the physics or 
chemistry in M33 may be different from that of AGNs and 
starbursts. The dominance of strong X-ray radiation in the nu- 
clei of AGNs or even of starbursts may be important for the 
differences in the line ratios, since it creates X-ray-dominated 
regions (XDRs) that change the chemical abundances (e.g. 
Meijerink & Spaans 2005i). 

The subsolar metallicity of M 33 may also play a role in cre- 
ating such a high HCN/HNC ratio. However, the HCN/HNC 
ratio obtained in M33 is significantly higher than those ob- 
served in similar low-metallicity environments, such as N159, 
30Dor in the LMC, as well as LIRS36 in the SMC, which are 
no hig her than 3.6 (besides a lower limit of 4.7 in N27 in the 
SMC) dChin et alJl998l:lHeikkila et al.lll999h and comparable t o 



the values found in Galactic GMCs ( Huettemeister et anil995h . 
Therefore, a subsolar metallicity alone does not seem to be a 
guarantee for very high HCN/HNC ratios. 

4.4.3. HCN vs. total infrared luminosity (Ltir) 

The Ltir/Lhcn ratios observed in M 33 (Table|2l) range between 
1300 and 3500 Lo/Kkms 'pc-. These ratios lie at the uppe r 
end of the values found in Milky Way clouds dWu et al.ll2010h . 
Normal galaxies show o n average total Lttr/Lhcn ratios of 
~ 900+70 Lo/Kkms ' pc^ dGracia-Carpio et al.l2008i GS04a,b). 
Higher values, in the ran ge of ~ 1100- 1700 Lp/K kms ' pc~, 
are found for (U-)LIRGs dGrac ia-Carpio et alJl2008l GS04a,b), 
while the highest reported ratios are found at the extreme end 
of the LIRG/ULIRG distribution, as well as in high-z galaxies . 



and range up to ~ 3900 Lp/Kkms ' pc^ dSolomon et a 
Gao et alj 120071; iGracia-Carpio et all 120081: IWuetaT 



2003 



201C; 



Garcfa-Burillo et akf 20121) . Thus the Lttp/Lhpm ratios in M33 
are among the highest ratios observed. 

In Sect.|6]below, we attempt to shed more light on the weak 
HCN emission and investigate the line ratios of CO, HCN, and 
HCO^ using models of photon-dominated regions that take not 
only the chemical network into account, but also the detailed 
heating and cooling processes of a cloud as well as radiative 
transfer. 

4.4.4. i-CO/'^^CO line ratio 

In our sample of GMCs in M 33 we find '^CO/'^CO line inten- 
sity ratios between 9 and 15. There is no obvious correlation 
with galacto-centric distance, FUV strengt h, or SFR. In a study 
of eight GMCs in the outer disk of M 33, iBraine et all d20ld) 
found similar ratios between 8.9 and 15.7. In the Milky Way, the 
isotopic ratio of '"C/'^C varies between values of 80 - 90 in th e 
solar neighborho od and 20 in the Galactic center dWilsonll 19991) . 
iPolk et all dl988h studied the '^CO/'^CO fine ratio in the Milky 
Way, for several large regions of the plane in comparison with 
the average emission from GMCs and of the centers of GMCs. 
They find that the average value rises from three in the centers of 
GMCs, to 4.5 averaged over GMCs, to 6.7 for the Galactic plane, 
with peaks of ~ 15. Their interpretation is that the higher ratios 
observed in the plane are caused by diffuse gas of only moder- 
ate optical thickness in '^CO. A similar interpretation may hold 
in M 33. The fraction of dense gas within the beam and optical 
depth effects may affect the ratios observed in M 33. 
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Table 4: Molecular abundances in M 33 and typical examples of 
galactic and extra-galactic sources. 



Source 


HCO+ 


HCN 


HNC 


C2H 


References 


M 33 min 


-10.3 


-9.8 


- 


- 




M 33 stackea 


-9.8 


-9.2 


-10.5 


-8.4 




M 33 max 


-9.5 


-9.1 








LMC N159 


-9.7 


-9.7 


-10.2 




1,2 


NGC253 


-8.8 


-8.3 


-9.0 


-7.7 


1,3 


M82 


-8.4 


-8.4 


-8.8 


-7.6 


1,3 


IC 342 GMC-A 






-8.7 


-7.5 


5 


Orion Bar 


-8.5 


-8.3 


-9.0 


-8.7 


1 


TMC-1 


-8.4 


-7.7 


-7.7 


-7.1 


1,6 


Transl. CI. 


-8.7 


-7.4 


-8.6 




1,4 



Notes. Entries show logiN{X)/N(H2)). 

Ref erences. (1) lOmontI ilOOH): ( 2) iJohansson et a ll (I1994D; 
(3) lMartmetal.1 ( l2006h : (4) iTumej (l2000h: (5) iMeier & Turner 



( 120051) and references therein; (6') lHogerheiide & SandeuTcOOQ) 



Ratios in the Magellanic clouds are observed by 
iHeikkila et all (1 19991) to cover values between 5 and 18, a 
somewhat wider range than found in M33. Unlike in M33, a 
gradient is seen on larger scales in a set of IR-bright nearby 
galaxies, dropping from values of about ten in th e center to 
values as low as t wo at larger radii dTan et al.ll201 ih . Although 
lAalto et aH (Il995h finds variations with galacto-centric radius 
in some galaxies of their IR-bright sample, other galaxies exist 
where the '^CO/'^CO stays constant with radius. They report 
a mean value of ~ 12 for the centers of most galaxies in their 
sample e xcept for the most luminous mergers with ratios of ^ 20 
(see also lCasoli et al.lll992h . 

5. Molecular abundances 

We use the observed line intensities of HCO^, HCN, HNC, and 
C2H to estimate the molecular column densities and abundances, 
assuming local thermodynamic equilibrium (LTE) and optically 
thin emission. Details on the calculations are given AppendixiDl 
with results shown in Table ID. fl The abundances are only lower 
limits in case the emission of HCN and HCO^ is optically thick. 
In comparison with our PDR-model analysis below we find, 
however, that for the best-fitting models to our observations the 
optical depths (t) in the centers of the lines of HCN and HCO^ 
are t < 0. 1, suggesting that emission is likely to be optically thin 
(cf. Sect. IT). In TablelH we compare the abundances with those 
found in other galaxies (LMC, NGC 253, M82, IC 342), in se- 
lected Galactic sources, the photon-dominated region Orion Bar, 
the dark cloud TMC-1, and a translucent cloud. The estimated 
column densities for the stacked values are given in Table |3] 

The abundances of HCO+, HCN, and HNC found in M33 
are very similar to those found in the LMC cloud N159. The 
abundances derived from the stacked spectrum of M 33 agree to 
within 0.5 dex with those of N159. Galactic sources have val- 
ues that are higher by more than an order of magnitude. The 
Orion Bar, for example, shows 1.8 dex to 0.8 dex higher abun- 
dances. The good agreement with the LMC may be driven by 
its s imilar metallicity o f 0.3-0.5 relative to the solar metallic- 
ity dHunter et al.ll2007|). which is onl y slig htly lower on aver- 
age than in M 33 jMagrini et a"i]|2007ll20Ta) . However, the C2H 
abundance observed in M33 and in the Orion Bar agree within 
0.3 dex. 

The LTE HCO+/HCN ab undance ratios in M 33 range be- 
tween 0.2 and 0.5 (Table lDTt . lGodard et alJ (l2010l) measured an 



HCO+/HCN abundance ratio of 0.5 in the diffuse ISM of the 
Milky Way, similar to the ratio found in the solar neighborhood, 
and similar to the higher ratios found in M 33. 



6. PDR models 

6.1. Setup 

To improve on the LTE analysis and to better understand why 
HCN is less luminous than HCO^ in M 33, we compare the ob- 
served HCO+/HCN, HCN/CO, and HCO+/CO line ratios with 
models of p hoton-d ominated regions (PD Rs) using th e Meudon 
PDR code jLe Petit et al. .2006; Gonzalez Garcia et al. ,20081) . 
The line intensities of the molecules '^CO, HNC, and C2H are 
not modeled. 

We ran a grid of models for different densities hh = 0. 1, 0.5, 
1, 5, 10, 50, 10- X 10'*cm-\ FUV fields Go = 10,50, 100 
in Habing unitfl and optical extinctions Ay - 2 -50 mag, i.e. 
in steps of log Ay ~0.2. We calculated this grid of models for a 
solar and a subsolar metallicity. The subsolar one is tailored to 
describe the metallicity in the disk of M 33. See Appendix lAlfor 
a detailed description of the model setup. 

6.2. Results 

The modeled HCN/CO Une ratio (Fig.|7] TablesEU and |R2] ) 
hardly varies with metallicity, «h or FUV. It varies, however 
strongly, with Ay. After an initial drop of upto an order of mag- 
nitude for extinctions less than about 4 mag, it rises until Ay 
~20 mag and then saturates at a value that is nearly independent 
of any of the input parameters. 

In contrast, the modeled HCO^/HCN ratio shows differences 
between the two metallicities for Ay <10mag. The subsolar 
model shows higher ratios than the solar model for a given Ay, 
FUV field, and density. Strong FUV fields and low densities in- 
crease the ratio. At higher optical extinctions the ratios are hardly 
influenced by changes in metallicity. Ay or FUV. This reflects 
the variation in HCN and HCO^ abundances. At low Ay, subso- 
lar HCN abundances are lower than solar ones by up to a fac- 
tor of 0.6 dex, while HCO^ is enhanced in the subsolar models 
by up to 0.7 dex. Also the CO abundances show a clear depen- 
dence on metallicity. For all input parameters, its abundance in 
the subsolar models is ~0.6dex lower than in the solar models. 
This directly reflects the underabundance of carbon of 0.6 dex in 
the subsolar models. As a result, the HCN/CO line ratio is fairly 
independent of metallicity. 

For optical extinctions in the range of Ay < 8 niag, w here 
the bulk of molecular gas in galaxies resides (lTielensll2005h . the 
HCN/CO ratio increases with increasing densities. In general in 
LIRGs/ULIRGs where most of gas has higher density than in 
normal galaxies, the HCN/CO is also higher (GS04a,b). 

7. Comparison with PDR models 

FigureQ shows the range in observed intensity ratios 
(cf. Table[2]i, where the plot shows that the optical extinc- 
tions play a decisive role in determining the line ratios. High 
extinctions of Ay > 16 mag are inconsistent with the observed 
HCN/CO ratios. Similarly, high densities of 10^ cm""* cannot 



- Habing units correspond to an average interstellar radiati on field 
(ISRF) between 6eV < hv < 13.6 eV of 1.6 lO^'ergcm-^ s"' dHabind 
196^). Another unit that is frequently used is the Draine unit of the local 
average ISRF, which is 1.7 x Go in Habing units. 
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Fig. 7: PDR model line ratios for subsolar (solid lines) and solar metallicity (dashed lines); HCO^/HCN (top) and HCN/CO (bottom). 
Different colors indicate different FUV field strengths Gq = 10 (red), 50 (green), and 100 (blue). Every panel of a subfigure shows 
the results for one density; from left to right and top to bottom rtu - 0.1, 0.5, 1, 5, 10, 50, and 10^ x lO'' cm"^. Gray areas mark the 
range of observed ratios in M 33. The dashed horizontal lines show the values from the stacked spectra. 



reproduce the measured HCO^/HCN ratios. Interestingly, 
metallicity only plays a minor role. In general, both metallicity 
models allow the observed range of line ratios to be reproduced. 

To quantify the agreement between the line ratios of the dif- 
ferent models and the individual observed clouds including the 
stacked spectrum, we use a x~ fit routine. To get a handle onto 
the errors of the best-fitting models, a Monte Carlo analysis is 
employed. The details on the fitting and Monte Carlo methods 
are given in Appendix |B] Table |5] shows the input parameters 
Ay, «H, and FUV of the best-fitting subsolar and solar metallic- 
ity models, i.e. those having the lowest for the stacked values 
and each observed cloud. 

7.1. Stacked ratios 

The best-fitting models for reproducing the stacked HCO^/HCN 
and HCN/CO ratios of 1.4 and 1%, respectively, that describe 
the averaged GMC properties are = 8mag, «h - 3 lO'^cm"-', 
and FUV = 68 Go . Emission stems from moderately dense gas 



with average line-of-sight column densities of 8 mag. The beam- 
filling factor Opuv derived from ratio of the beam-averaged TIR 
intensity to the fitted local FUV field is -30%; i.e., the fitted 
FUV field strengths are significantly higher than expected from 
the observations. The same holds for the beam filling factors de- 
duced from the ratios of extinctions derived from CO and the Ay 
of the best-fitting models, which are about O^,, = 50%. This is 
not surprising, however, and indicates that emission is clumped 
within the 1 14 pc beam. From the models of the calculated grid 
closest to the best fit, i.e., Av = 6 and 10 mag, «h = 1 10"*, 
and FUV = 50 Go, we find optical depths (t) in the centers of 
the lines of HCN and HCO+ that He between 0.02-0.1 and 
0.07 - 0.1, respectively. Thus both lines are optically thin. '^CO 
is moderately optically thick with optical depths of t 4 - 25. The 
line width assumed in the Meudon PDR code is ~ 3kms ' (cf. 
AppendixlAli. 
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Table 5: Best-fitting PDR models 





HCO+/HCN 


HCN/CO 




"H 


FUV 




a 












[%] 


[mag.] 


[cm-^] 


[Go 


] 














subsolar metallicity Models 


Stacked 


1.4 ± 


0.2 


1.0 ±0.1 


8±3 


(3 ± 4) 10" 


68 ± 


24 


0.5 ± 


0.2 


0.3 


±0.1 


0.1 


±0.1 


N06 


2.5 ± 


0.2 


1.4 ±0.3 


4±2 


(6 ± 4) 10^ 


41 ± 


16 


1.7 ± 


1.1 


0.9 


± 0.4 


1.4 


±0.5 


N03 


1.9 ± 


0.1 


0.8 ± 0.2 


8 ± 2 


(3 ± 2) 10' 


27 ± 


20 


0.7 ± 


0.2 


0.8 


±0.6 


0.6 


±0.3 


GMC91 


1.4 ± 


0.3 


0.4 ±0.1 


6±2 


(1 ± 5) 10" 


54 ± 


44 


0.9 ± 


0.3 


0.2 


± 0.2 


0.1 


±0.1 


NO! 


1.4 ± 


0.2 


1.7 ± 0.4 


7±4 


(3 ± 4) 10" 


42 ± 


40 


0.3 ± 


0.2 


0.3 


±0.3 


0.2 


±0.1 


GMCl 


1.1 ± 


0.1 


2.9 ± 0.4 


10 ± 3 


(1 ± 2) 10" 


30 ± 


37 


0.6 ± 


0.2 


1.7 


± 2.1 


0.7 


±0.5 


GMC26 


1.2 ± 


0.1 


1.0 ± 0.2 


9±2 


(1 ± 3) 10" 


67 ± 


24 


0.3 ± 


0.1 


0.2 


±0.1 


0.1 


±0.1 


solar metallicity models 


Stacked 


1.4 ± 


0.2 


1.0 ±0.1 


10 ± 1 


(5 ± 1) 10' 


47 ± 


12 


0.4 ± 


0.0 


0.5 


±0.1 


0.6 


± 0.4 


N06 


2.5 ± 


0.2 


1.4 ±0.3 


2± 


(5 ± 1) 10' 


100: 


t4 


3.1 ± 


0.8 


0.4 


±0.0 


3.4 


± 1.2 


N03 


1.9 ± 


0.1 


0.8 ± 0.2 


6± 3 


(2 ± 2) 10' 


81 ± 


32 


0.9 ± 


0.5 


0.3 


±0.1 


3.4 


±0.9 


GMC91 


1.4 ± 


0.3 


0.4 ±0.1 


4± 1 


(1 ± 1) 10' 


98 ± 


11 


1.4 ± 


0.3 


0.1 


±0.0 


0.5 


± 0.4 


NOl 


1.4 ± 


0.2 


1.7 ± 0.4 


7±4 


(5 ± 3) 10' 


50 ± 


43 


0.4 ± 


0.2 


0.3 


± 0.2 


0.2 


± 0.2 


GMCl 


1.1 ± 


0.1 


2.9 ± 0.4 


11 ± 2 


(3 ± 6) 10' 


26 ± 


34 


0.6 ± 


0.1 


1.9 


± 2.6 


0.5 


± 0.2 


GMC26 


1.2 ± 


0.1 


1.0 ± 0.2 


9±2 


(5 ± 1) 10' 


52 ± 


17 


0.2 ± 


0.1 


0.2 


±0.1 


0.4 


±0.3 



Notes. beam-filling factor derived from Av(''CO)/Av(Model); beam-filling factor derived from FUV(TIR)/FUV(Model); average of the 
best-fitting models. 



7.2. Individual regions 

Here, we focus on individual regions grouped by their particular 
HCO^/HCN ratios and thus their best-fit Ay values: no6 show- 
ing a high ratio of 2.5, GMC91, no3, and nol have intermediate 
ratios of 1.4-1.9 and GMC26, GMCl having ratios of 1.1-1.2. 

GMCno6 This cloud shows the highest HCOVHCN ratio of 
2.5, while at the same time the HCN/CO ratio is relatively low 
with 1 .4 and weaker than expected from the linear fit to the M 33 
data (cf. Sect. l4.4.1l and Fig.©. GMC no6 is best fitted by subso- 
lar models that yield a low best-fitting value for Ay of ~ 4 mag, 
while the best-fitting density and FUV strength are 6 10^ cm""* 
and 40 Go, respectively. The beam-filling factor derived from Ay 
is 1.7, indicating that emission completely fills the beam with 
several clouds along the line-of-sight. This cloud has the second 
highest star formation rate of our sample of 35.9 Mq Gyr"' pc ~, 
and the same holds for the FUV field strength of 37.3 Gq. 

GMC91, no3, and nol The line ratios of these three clouds 
are best described by subsolar models. The best-fitting Ay are 
similar with 6-8 mag. So are the FUV 30-50 Go and the den- 
sities 3 lO^'cm"-' - 3 10"* cm"^*, and nol and no3 have similar 
SFR rates of ~ 13 M0 Gyr"' pc"^, while no3 is a factor of four 
more massive than nol with Mh, =810^ Mq. GMC91 lies at 
only 320 pc distance in close vicinity of GMCno3 and is only 
slightly more massive than the same. It is the most CO intense 
cloud in our sample while its HCN and HCO^ emission is rel- 
atively weak. This renders GMC91 somewhat peculiar and re- 
sults in a low HCO^/CO ratio of 0.6% and, as already found 
by RPGl 1, a particularly low HCN/CO intensity ratio of 0.4 % 
(Table|2|. The HCN/CO ratio of GMC91 is much lower than the 
ratios observed in the disk of the Milky Way bv iHelfer & BUtd 
(Il997h . who find 2.6 % ± 0.8 % and also in the inner disks (5- 
lOkpc) of normal galaxies by GS04b who find 4 % ± 2 %. The 
relatively weak HCN and HCO^ emission may indicate a low 
fraction of dense mass in GMC91, which thus may be a rather 
quiescent GMC with only a low SFR. Indeed, it has the lowest 
star formation rate of 4 M© Gyr"' pc"^ of all observed clouds. 



GMCl and GMC26 These two clouds host the lowest observed 
ratios of HCO^/HCN of 1.1 - 1.2. Here, the solar models provide 
slightly better or equal fits than the subsolar models. However, 
since the ISM of M 33 is subsolar, here we discuss only the best 
fits to the subsolar models. These GMCs have similar best-fitting 
input parameters of Ay 9-10 mag and «h ~ lO'^cm"^, while the 
best-fitting FUV field strengths are ~ 30 Go and 70 Go, respec- 
tively. However comparing their physical properties in Table |2] 
again, these two clouds are actually not at all alike. GMCl is 
located in the very center of M 33 and is by far the most mas- 
sive cloud in our sample. It is actively forming stars at a rate of 
65 M© Gyr"' pc"^, the highest in our sample, and has an corre- 
spondingly high FUV field of 50.7 Gq. GMC26 has much lower 
HCN/CO ratios, and its SFR is only 6.6 M© Gyr ' pc"^ the sec- 
ond lowest in the sample with exception of no2. 

For all best-fitting solutions of the six individual positions 
we find that the modeled optical depths of HCN and HCO^ are 
T < 0.1, which renders emission of these lines to be likely opti- 
cally thin. This also justifies the assumption of optical thin emis- 
sion for the UTE analysis in Sect|5] Indeed, the PDR modeled 
abundances of both molecules are comparable to the ones de- 
rived from LTE (cf. Table|4]and|E2]i. 

In conclusion, it is noteworthy to repeat that the line ratios 
studied here are fairly independent of the metallicity, SF activity, 
and FUV field strength of the parent GMC, while the optical 
extinction has a major influence on the modeled line ratios. 

8. Summary 

We present IRAM 30m observations of the ground-state tran- 
sitions of HCN, HCO+, '^CO and '^CO of seven GMCs dis- 
tributed along the major axis in the disk of the nearby spiral 
galaxy M33. We achieve a spatial resolution of ~ 114pc at a 
frequency of 89 GHz. 

The molecular gas masses of the target GMCs vary by a 
factor of ~ 130 between 0. 1 10^ Mq (GMC no2) and 13 10^ Mq 
(GMC 1 ) and the star formation rates derived from Ha and 24 [a.m 
images vary by a factor of more than 50. The FUV field strengths 
show a variation of more than a factor 20. Below, we summarize 
the main results. 
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For the six GMCs where HCO^ is detected, peak hne tem- 
peratures (on the scale) vary between 6 and 12 mK. The 
HCO^/HCN-integrated intensity line ratios lie between 1 . 1 
and 2.5 (on the Kkms ' scale, cf. Table d. Similar line r a- 
tios are observed in the disk of M 3 1 GBrouillet et al.ll2005l) . 
The line intensity ratios HCN/CO and HCO+/CO vary 
between (0.4-2.9)% and (0.6-3.5)%, respectively. The 
spread of ratios foun d in M 33 is slightly l arger than in the 
spiral arms of M31 dBrouillet et alJl2005L Fig|5]). GMC91 
exhibits a particularly low HCN/CO ratio of 0.4 %, which is 
much lower than valu es in the Galactic disk of 2.5 % + 0.6 % 
dHelfer & Bhtz|[T997l) or in normal galaxies with 4 % + 2 % 
(GS04a). 

The LxiR/i'HCN luminosity ratios range between 1.3 10^ and 
3.5 10^ an d are situated at the very high end of ratios found 
bv IWu et' al. (2010) in molecular clouds of the Milky Way 
and LIRGs/ULIRGs. This shows that HCN emission in com- 
parison to the LxiR in M 33 particularly weak. 
Stacking of all spectra taken at the seven GMC positions 
leads to 3 cr detections of CCH and HNC. The HCN/HNC 
ratio of 5.8 is remark able high. It is higher than val- 
ues found in the LMC dChin et al.lll997lll998h . in IC 342 
([Meier & Turiieil I2005h. in samples of LIRGs/ULIRGs 



( Baan et al. 20081: Costagliola et al. 201 Ih . in starburst and 
Seyfert galaxies (e.g. 'Aalto et al.' '2002*), and in Galactic 
molecular complexes [Wootten et aLi (il978i) . where all to- 
gether no values higher than three are reported. 
The HCO^, HCN, HNC abundances, derived assuming LTE, 
agree with those of the LMC cloud N159 within 0.5 dex. In 
contrast, the Orion Bar, a Galactic massive star-forming re- 
gion, shows significantly higher abundances of all three trac- 
ers by 0.8 dex to 1.8 dex. These striking differences may re- 
flect the factor two subsolar metallicities of both the LMC 
and M 33. 

Employing the Meudon PDR code to model photon- 
dominated regions we investigated the influence of the 
metallicity on the abundances and emission of HCN and 
HCO^. For a range of optical extinctions, volume densi- 
ties, and FUV radiation field strengths, we derived two 
sets of models with different met allicity, one reflecting the 
abundances in the Orion nebula bv lSimon-Diaz & Stasiriskal 
( 201 lb. the other th e average subsolar metalUcity of M33 
dMaerini et alJl2010h . 

Both sets of models are able to describe the observed range 
of HCO+/HCN and HCN/CO line ratios reasonably well 
(X^ < 3.4). Therefore, changes in metallicity do not need 
to be invoked to describe the observed line ratios. The obser- 
vations are described by subsolar models with optical extinc- 
tions between 4 mag and 10 mag and moderate densities of 
<3 10"* cm"^, with little influence by the FUV field strength. 
The optical extinction has a pronounced influence on the 
modeled ratios, while FUV field, metallicity and even den- 
sity only play minor roles. The modeled lines of HCN and 
HCO^ of the best-fitting models are found to be optically 
thin with optical depths t < 0. 1 . 
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Table A. 1 : Initial Abundances used in the PDR models 



Species 


Sun" 


Orion 


M33^ 


diff. in [dex]'' 


TT 

til 


U.l 


U.l 


U.l 




H 


0.8 


0.8 


0.8 




He 


0.1 


0.1 


0.1 




C 


8.43 


8.37 


7.77 


0.60 


O 


8.69 


8.65 


8.27 


0.38 


N 


7.83 


7.92 


7.31 


0.61 


S 


7.12 


6.87 


6.71 


0.16 


Fe 


7.50 


6.00 


5.73 


0.27 



Notes . En tries show 12+ lo,g(;i(X)/( 2«(H2 )+«( H)). |A splund etal] 
(120091): ISimon-Di'az & Stasinska l l2011h : iMagrini et al.i COloF : 
iHenrv et al.l ( 120001) : between Orion and M 33. 

Appendix A: PDR Model setup 

We use the Meud o n PD R code (iLe Petit et alj 120061: 
iGonzalez Garcia et al. I l2008l) . which solves the thermal 
balance, chemical network, detailed balance, and radiative 
transfer of a plane-parallel slab of optical extinction Ay and 
constant density hh illuminated from both sides by an FUV 
field of intensity Go. The modeled line widths are calculated by 
the Meudon code via the Doppler broadening of the gas due to 
the kinetic temperature and the turbulent velocity of the gas. 
The latter dominates the line widths and is an input parameter 
to the Meudon code. We use the default value of the Meudon 
code which, expressed in terms of the FWHM of the velocity 
distribution, is (A v)fwhm ~ 3 km s"' . This is slightly lower than 
the observed line widths that are in the range of 4-11 kms"'. 
We also used the defaults for additional input parameters (e.g. 
the cosmic-ray ionization rate and size distribution of dust 
grains). 

The grid of models has been calculated for volume densities 
nil = 0.1, 0.5, 1, 5, 10, 50, 10^ x 10'*cm-^ FUV field strengths 
Go = 10,50,100 in Habing units, and optical extinctions 
Ay = 2 -50 mag, i.e. in steps of log Ay ~0.2. 

The input range of FUV fields is motivated by the range of 
FUV fields derived from the TIR in the six observed clouds 
which vary between Go = 11.3 and 50.7, excluding position 
no2 with a very low Go of 2.5 (cf. Table |2]i. We add mod- 
els of Go =100 to cover possible higher local radiation fields. 
The values of Ay range from those found in translucent clouds 
with 2 mag to extinctions of 50 mag typically found in resolved 
Galactic high-mass star-forming clouds. The densities cover val- 
ues found in a typical molecular cloud, covering the critical den- 
sities of the observed tracers and transitions. 

We calculated this grid of models for solar and subsolar 
metallicities. Table lA. ll shows the abundances measured in M 33 
and in the Orion nebula which we use as initial abundances for 
our PDR models. The abundances of the sun are listed for com- 
parison. For the "solar" metallicity model, we adopt the abun- 
dances of O, C, N, S, Fe measured in the H ii region of the Orion 
nebula bv lSimon-Dfaz & Stasinskal (1201 lb . The heavy elements 
S and Fe are depleted to dust grains with respect to the abun- 
dances found in the solar photosphere. 

For the subsolar metallicity PDR model, we adopt the 
averaged O, N, an d S abundances measured in M33 from 
IMagrini et al1(l2010l) who targeted 33 H ii regions between 1 and 
8 kpc galacto-centric distance. The carbon gas-phase abundance 
has not been measured in M 33. We estimate it using t he metal- 
licity dependence of the C/O ratio described by CHenrv et alJ 
I2OOOI) . To derive the Fe abundance, we scale the solar Fe/O ratio 
by the subsolar O metallicity of M 33. 



Appendix B: Chi-squared fitting and IVIonte Carlo 
analysis 

To quantitatively compare the modeled and observed line ratios 
we calculate the summed squared residuals weighted by the ob- 
servational errors between the modeled and observed line ratios 
HCO+/HCN and HCN/CO, for every model and for each ob- 
served cloud as 

2 1 (RatiOobs,i - RatiOmod,i)^ 



where cr is the observational error of the line ratios deduced from 
the errors of the line intensities, cf. Sect.|2l We normalize thex~ 
by dividing by = 2, the number of independent observed line 
ratios. The minimum gives the best-fitting model. 

To get a handle on the errors of the best-fitting models, we 
employ a Monte Carlo simulation. We assume that the measured 
ratios and their errors follow a Gaussian distribution with an ex- 
pected value equal to the ratio and a variance equal to cr^. 

We generate 5000 sets of HCN/CO and HCO+/HCN ratios 
drawn randomly from their Gaussian parent distribution. For ev- 
ery set we calculate ax^ as explained above and deduce the best 
fits. From all best fits of all sets we derive the mean and the 
standard deviation of the input parameters: nn. Ay, FUV, andx^ 
(Table|5]l. 



Appendix C: GMC properties 

For comparison we derived complementary properties, the star 
formation rate, the FUV field, the Xco factor, the masses of the 
molecular and atomic material, and the optical extinction within 
the beam at the observed positions (Table|2]i. Below, we explain 
in detail how these values have been derived. All errors are cal- 
culated using Gaussian error propagation. The particular uncer- 
tainties of individual measurements used to derive the comple- 
mentary quantities are given below. 

C.1. Star formation rate (SFR) and total infrared luminosity 
(TIR) 

To obtain an ex tinction free tracer of the SFR we com- 
bine the 24 \im ("Tabataba ei et al.l l2007h and Ha luminosity 
(iHoopes & Walterbos 200 0) using SFR = 5. 3 IQ-'^^lLiHoc) + 
(0.031 + 0.006)L(24 ^m)] fcalzetti et al.ll2007h . We assume er- 
rors of 7% for 24 [a.m (Spitzer Observers Manual v8.0) and 15% 
for Ha. 

Another measure of the SFR is the total infrared luminosity 
(Ltir) emitted between 1 fim and 1 mm. We estimate it by adding 
the measured brightnesses at 24, 70, 100, 160, 25 350, 500 urn 
in W k pc~^ weighted by the factors determined bv lBoquien et al.l 
(1201 11) . We convert the determined brightnesses to luminosities 
corresponding to the 28 " beam by multiplying with the beam 
size in kpc T he errors on the deriv ed TIR are the cr values of 
the fits given bv lBoquien et al.l (1201 ih . 

C.2. Molecular line luminosities 

We calculate the molecular line luminosities L[^^j^ and ij^co* 
units of Kkms"' pc^ by multiplying the integrated line intensi- 
ties with the beam size in pc^, i.e. 1 .5 10"^ pc^. Errors follow from 
the measurement uncertainties given in Table |2] 
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Fig.A.l: '-CO and CO 1 -0 spectra at the seven selected GMC positions at high-velocity resolution of 1 kms '.The '-CO spectra 
of positions GMC91, GMC26 and GMCl are taken from RPGll. The blue lines overplotted on the '^CO spectra show the Gaussian 
line fits used to deduce the FWHM of the lines which are listed in Tabled 



C.3. FUVfluxGo 

The FUV flux Go impinging on the dusty cloud surfaces can 
be estimated from the emitted total infrared intensity assum- 
ing that all FUV photons are absorbed by the d ust. We use 
Go = 47r0.5LTiR (cf. e.g. iMookeriea et all l20Hl) with G o in 



units of the Habing field 1 .6 10 ^ erg s ' 



cm (lHabingll968h and 

the TIR intensity in units of erg s ' cm"^ sr '. Errors follow di- 
rectly from the errors on Ltir- The FUV flux that we use to com- 
pare with values deduced from the stacked spectrum is the aver- 
age FUV flux of all individual clouds, i.e. 21.6 + 0.7Go. 



C.4. Xco factor and atomic and molecular gas mass 

To derive the masses of the molecular gas we first convert the 
CO 1-0 integrated intensities into H2 column densities (A^hO via 
A^Hi = ^co^co- We assume that the conversion factor Xco 
(c f. Table [2I1 is a function of total infrared luminosity as found 
bv llsraell(ll997h for a sample of low-metallicity envi ronments in 
the L MC and SMC and for two clouds in M 33 (cf. iLerov et all 
1201 Ih . 

[Israeli (1 19971) derived A^h, for CO clouds in the LMC, the 
SMC, and other low-metallicity galaxies from their far-infrared 
surface brightness and Hi column densities. He finds local Xco 
values up to a magnitude higher than the widely used Gala ctic 
value of 2 lO^o (K km s ') ' cm ^ (e.g. iDame et al.ir200lh . In 
particular he studies the two bright H 11 regions NGC 604 and 
NGC 595 in M 33, for which he derives X-factors of 22 10^° (K 



km s-')-' cm-^ and 12 10^" (K km s"') 



l llsraell(ll997h 



further finds a correlation between Xco and the TIR-flux. For 
NGC 604 we find a TIR luminosity of 27.18 and 7.82 10*= L© 
for NGC 595 (cf. Sec lCTt . We estimate the local Xco for our 
seven clouds from their TIR by using the mean of the Xco/i^TiR 
ratio of NGC 604 and NGC 595 as a conversion factor, i.e. Xco 
= 1.17 10-" 10 Lq Ltir- The derived X-factors vary between 7 
and 1.5 in units of 10^° (K km s ') ' cm~^ for the six positions 
with HCN and HCO+ detections (Table|2ll. iGratier et al.] (1201 Ol) 



,-2 



assume a constant X-factor of 4 10 (K km s ) cm 

The mass of the molecular gas Mh, is estimated from the in- 
tensity of the '^CO 1-0 line using the Xco-factor derived in the 
previous section and Mh, - 2x 1.36 A^h, nip A with the proton 
mass mp, the beam size A. The factor two accounts for the two 
protons of molecular hydrogen and 1.36 includes the contribu- 



tion of heavy elements to the total mass. Errors follow from the 
measurement uncertainty of '^CO given in Table |2] 

The mass of the atomic gas is estimated from the intensity 
of die Hi Une via the column density A^hi = 1-82 x 10"^ /hi 
(iRohlfs & Wilsonll2000l) assuming optically thin emission. The 
atomic gas mass is calculated with the same formula as for the 
molecular mass without the factor two accounting for the two 
protons in H2. Errors assume a 15% calibration error on the HI 
data. 



C.5. Optical extinction Ay 



As shown by iBohlin et al.l (Il978l) a correlation between the 
amount of hydrogen atom column density N(Hi + H2) and 
the color excess exists; they estimate a conversion factor of 
<N(Hn-H2)>/E(B-V)= 5.8 10^' atoms cm'^ mag' '. The color 
excess and the optical extinction Ay are also linked. iBohlin et all 
(Il978h give a value for the interstellar average of R = Av/E(B- 
V) = 3.1. Using the two conversion formula and the values for 
Nhi and Nh,, derived as explained in the previous section, we 
calculate an estimate of the optical extinction for the observed 
clouds per beam. For clouds not filling the 115pc beam, local 
extinctions may be higher. The errors given for Ay (cf. Table |2]i 
include the errors of N(Hi) and N(H2). We estimate the Ay for 
the stacked spectra by averaging the values of the individual po- 
sitions, i.e. 4.2 + 1.9 mag. 



Appendix D: LTE column densities 

We estimate column densities from the integrated intensities, as- 
suming UTE and optically thin emission: 



N : 



3/1 Z expiw;;) 



87rV- J [1 



■ expi-^)] 



{Jv(Tex)-Jv(TBc)) 



with: 



liv 



1 



T„,,,dv,(D.l) 



(D.2) 



where Z the partition function, J the rotational quantum num- 
ber of the upper level, /i the dipole moment, Tex the exci- 
tation temperature, Tbg the temperature of the cosmic back- 
ground, and the main beam temperature of the specific line. 
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Table D. 1 : LTE column densities at the individual positions. 
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\7fU r<-»m~^l'\ T TP 
/V^^rl2 [cm \) - L,lil 


1 IQaiOl 

1 . loe+zi 






1 Q/1o_i_01 


4. 14e+Zi 


V.ZoC+ZU 


j.Zie+ZU 








Q ^9p4-9n 


9 ASp4-9 1 


9 SSp-i-91 


U.J ICT^ZU 


zL nnp-4-1 


log yV(''CO(l-0))/A'(H2) 


-5.00 


-5.22 


-5.07 


-5.09 


-5.07 


-5.14 


-5.15 


log 7V("CO(2-l))/A'(H2) 


-5.31 


-5.48 


-5.46 


-5.50 


-5.51 


-5.36 


-5.80 


log A'("CO(1-0))/A'(H2) 


-5.85 


-5.85 


-5.85 


-5.85 


-5.85 


-5.85 


-5.85 


log A'(HCO+)/A'(H,) 


-9.46 


-9.74 


-10.01 


-9.91 


-10.34 


-9.78 


< -10.15 


log iV(HCN)/yV(H2) 


-9.14 


-9.07 


-9.37 


-9.48 


-9.79 


-9.20 


-9.09 


iVCHCO+VyVCHCN) 


0.48 


0.21 


0.23 


0.37 


0.28 


0.26 


<0.09 


Tex [K] 


23 


25 


22 


22 


20 


23 


19 



Table D.2: Molecule parameters 





y 










[GHz] 


[Debye] 


[K] 


[cm-3] 


C2H(l-0, 3/2-1/2) 


87.2841050 


0.77 


4.19 


1 10' 


HCN 1-0 


88.6316022 


2.99 


4.25 


2 10^ 


HCO+ 1-0 


89.1885247 


3.90 


4.28 


3 10" 


HOC+ 1-0 


89.4874140 


2.77 


4.29 


3 10" 


HNC 1-0 


90.6635680 


3.05 


4.35 


2 10^ 


"CO 1-0 


110.2013543 


0.11 


5.28 


4 10^ 


i^CO 1-0 


115.2712018 


0.11 


5.53 


4 10^ 



Notes. for collisions with Hi neglecting opacity effects. 



Table [De lists frequency, fiu/k*, and the critical densities of 
the observed transitions. The partition function is taken from 
the "Cologne Database for Molecular Spectroscopy" (CDMS) 
jMiilleret al.li2001, .2005i) . The value of Z depends on Tex. The 
CDMS lists values for discrete temperatures in the range from 
2.75 to 500 K for every molecule. We use these values to inter- 
polate Z for the excitation temperature we adopt at a specific 
position. 

To estimate the excitation temperature, we derive the dust 
temperature, assuming that Tex, the kinetic temperature Tkin, 
and the dust temperature Tdust are coupled and equal. To 
estimate Tdust we obtain the dust SEDs for all our point- 
ings from the Herschel and Spitzer observations at 500, 
350, 250, 160, 100, and 24 yum 28"resolution and fit two- 
component graybody models to the data. Each component is 
described by S'v = B{v, T)Ty = B{v, T)KyMd / D~ , assuming op- 
tically thin emission, with the flux S,, , the Planck function 
By , the opacity Ty , the dust mass , the distance D, and 
the dust absorption coefficient K y - .4(v/(250 GHz)/ cm^ g"' 
dKruepel & Siebenmorgerilll994t iKrligel 2003]), /3 is the dust 
emissivity index assumed to be 1.5. 

Tdust is found to vary only slightly, between 20 K in the 
northernmost position observed, GMC no2, and 26 K in the nu- 
clear GMC, GMCl (cf. Table IdTTTi. We assume that Tex equals 
the temperature of the cold dust component. 

The CCH 1-0 transition is split into six hyperfine compo- 
nents. In the stacked spectrum, we detected the 3/2-1/2, 2-1 
transition, which is the strongest one. To determine LTE column 
densities of CCH, we estimate the total integrated intensity of 
the CCH 1-0 line including all hyperfine components by divid- 
ing the result of a Gaussian fit to the detected co mponent by its 
relative strength, i.e. 0.416 dPadovani et alJl2009l) . 



The column density of H2 is estimated from CO 1-0 
, assuming a '~CO/'^C O abundance ratio of 60 as typically 
found in the Milky Way (iLanger & Penziasll 19931) and a H2/CO 
abundan ce ratio of 8.5 10"^, also found for the Milky Way 
(Frerking et al.l [19821) . From this we derive the relative abun- 
dances of the different species with respect to A^h, ■ The derived 
column densities, as well as the relative abundances, are given 
in Table ID. ll for individual clouds and in Table[3]for the stacked 
spectra. For comparison we give in Table ID. li the values of A^h^ 
that we derive from the '^CO 1-0 lines using the individual Xqo- 
factors derived for every cloud (see above). After comparing 
both methods to obtain the column density of H2, we find that 
the results are consistent within a factor of 2.5 excluding posi- 
tion no2. The values of /Vh, for the stacked values in Table[3]are 
deduced using an average over the individual Xco factor of the 
observed GMCs of 3.0 10^° (K km s^i)-' cm"". 

Appendix E: PDR Model results 

The modeled abundances and intensity ratios of the ground- 
state transition of HCN, HCO, and '^CO using the Meudon 
PDR code for subsolar and solar metallicities are shown in 
Tables IE. II and IE.2I Both tables are subdivided in three blocks 
one for each modeled radiation field of Go = 10, 50, and 100 
in Habing units, which illuminates the modeled clouds from 
both sides. The other two input parameters that have been var- 
ied are the density with values of hh - 0.1, 0.5, 1, 5, 10, 
50, 10" X lO'^cm"^ and the optical extinctions with values of 
Ay = 2, 4, 6, 10, 16, 26, 40, and 50 mag. 
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Table E.l. Results of the PDR models with subsolar metallicities for M 33 (cf. Sect.[6ll. 
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/I A A 

4U.U 


^ in3 
3 10 


O A 


-8.4 


A A 

-4.U 


11 A A 

1 1.44 


1 O AO 

iz.yo 


1.14 


/I A A 

4U.U 


1 1 n4 
1 10 


A 1 

-y.V 


O A 

-8.4 


/I A 

-4.U 


1 1 CA 

1 i.jy 


11 11 

11.11 


A AA 

u.yo 


40 






O.J 


^ Q 


1 j.Uj 


f\ 7Q 

O. / 7 




40.0 


1105 


-9.8 


-8.6 


-3.9 


14.05 


5.25 


0.37 


40.0 


5105 


-10.2 


-9.0 


-3.9 


10.59 


2.63 


0.25 


40.0 


1 10* 


-10.3 


-9.2 


-3.9 


8.08 


1.86 


0.23 


50.0 


1 103 


-8.2 


-8.1 


-4.0 


15.51 


19.81 


1.28 


50.0 


5 103 


-8.8 


-8.3 


-4.0 


11.89 


14.17 


1.19 


50.0 


1 10* 


-9.1 


-8.3 


-4.0 


12.23 


12.53 


1.02 


50.0 


510* 


-9.6 


-8.4 


-3.9 


14.97 


8.28 


0.55 


50.0 


1105 


-9.8 


-8.6 


-3.9 


15.82 


6.54 


0.41 


50.0 


5105 


-10.1 


-9.0 


-3.9 


12.80 


3.35 


0.26 


50.0 


1 10* 


-10.3 


-9.2 


-3.9 


9.96 


2.39 


0.24 


Go = 50 


2.0 


1103 


-11.1 


-13.0 


-8.7 


1.19 


219.55 


185.29 


2.0 


5103 


-11.3 


-11.7 


-7.6 


2.38 


15.95 


6.71 


2.0 


1 10* 


-11.4 


-11.2 


-7.1 


2.29 


4.48 


1.95 


2.0 


5 10* 


-11.4 


-10.4 


-5.8 


0.98 


0.37 


0.38 


2.0 


1105 


-10.9 


-10.2 


-5.2 


0.53 


0.33 


0.61 


2.0 


5 105 


-11.0 


-9.6 


-4.5 


0.27 


0.04 


0.13 
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Table E.l. Continued 







A'(HCO+)° 


X(HCN)° 


XCCO)" 


HCN/CO 


HCO+/CO 


HCO+/HCN 


(mag.) 


(cm-3) 








% 


% 




2.0 


1 10" 


-11.1 


-9.4 


-4.3 


0.29 


0.02 


0.06 


4.0 


1103 


-11.0 


-11.6 


-6.7 


0.24 


2.34 


9.64 


4.0 


5103 


-11.2 


-10.8 


-5.3 


0.11 


0.09 


0.75 


4.0 


1 10" 


-11.3 


-10.5 


-4.9 


0.10 


0.04 


0.40 


4.0 


510" 


-11.2 


-10.0 


-4.4 


0.18 


0.03 


0.16 


4.0 


1105 


-11.0 


-9.8 


-4.3 


0.24 


0.04 


0.18 


4.0 


5105 


-11.2 


-9.5 


-4.1 


0.40 


0.02 


0.05 


4.0 


1 10" 


-11.3 


-9.4 


-4.1 


0.45 


0.01 


0.03 


6.0 


1103 


-10.8 


-11.0 


-5.8 


0.20 


0.57 


2.81 


6.0 


5103 


-11.0 


-10.4 


-4.8 


0.14 


0.09 


0.59 


6.0 


1 10" 


-11.0 


-10.2 


-4.5 


0.18 


0.07 


0.38 


6.0 


510" 


-11.1 


-9.8 


-4.2 


0.34 


0.05 


0.16 


6.0 


1105 


-11.0 


-9.7 


-4.1 


0.41 


0.06 


0.15 


6.0 


5105 


-11.2 


-9.5 


-4.1 


0.57 


0.03 


0.05 


6.0 


1 10" 


-11.3 


-9.4 


-4.0 


0.59 


0.02 


0.03 


10.0 


1103 


-10.2 


-10.0 


-5.1 


1.06 


1.41 


1.33 


10.0 


5103 


-9.8 


-9.6 


-4.3 


1.02 


1.54 


1.51 


10.0 


1 10" 


-9.9 


-9.7 


-4.2 


0.84 


1.08 


1.29 


10.0 


510" 


-10.3 


-9.7 


-4.1 


0.74 


0.48 


0.65 


10.0 


1105 


-10.4 


-9.6 


-4.0 


0.75 


0.36 


0.48 


10.0 


5 105 


-10.8 


-9.6 


-4.0 


0.77 


0.14 


0.18 


10.0 


1 10" 


-10.9 


-9.5 


-4.0 


0.74 


0.09 


0.12 


16.0 


1103 


-8.7 


-8.4 


-4.3 


13.57 


14.65 


1.08 


16.0 


5103 


-9.2 


-8.8 


-4.1 


6.49 


6.21 


0.96 


16.0 


1 10" 


-9.4 


-8.8 


-4.1 


5.54 


4.16 


0.75 


16.0 


5 10" 


-9.9 


-9.0 


-4.0 


4.23 


1.80 


0.43 


16.0 


1105 


-10.0 


-9.1 


-4.0 


3.43 


1.31 


0.38 


16.0 


5105 


-10.4 


-9.4 


-4.0 


1.79 


0.55 


0.31 


16.0 


1 10" 


-10.5 


-9.5 


-4.0 


1.38 


0.36 


0.26 


26.0 


1103 


-8.4 


-8.2 


-4.1 


15.23 


17.30 


1.14 


26.0 


5103 


-9.0 


-8.5 


-4.0 


9.56 


9.90 


1.04 


26.0 


1 10" 


-9.2 


-8.5 


-4.0 


8.88 


7.56 


0.85 


26.0 


5 10" 


-9.7 


-8.7 


-4.0 


8.42 


3.97 


0.47 


26.0 


1105 


-9.8 


-8.8 


-4.0 


7.58 


2.97 


0.39 


26.0 


5105 


-10.2 


-9.2 


-4.0 


4.22 


1.30 


0.31 


26.0 


1 10" 


-10.3 


-9.3 


-4.0 


3.02 


0.85 


0.28 


40.0 


1103 


-8.3 


-8.1 


-4.0 


15.47 


18.89 


1.22 


40.0 


5103 


-8.9 


-8.4 


-4.0 


10.74 


12.44 


1.16 


40.0 


1 10" 


-9.1 


-8.4 


-4.0 


10.36 


10.40 


1.00 


40.0 


5 10" 


-9.6 


-8.6 


-4.0 


11.18 


6.47 


0.58 


40.0 


1105 


-9.8 


-8.7 


-4.0 


10.91 


5.01 


0.46 


40.0 


5105 


-10.1 


-9.1 


-3.9 


6.97 


2.27 


0.33 


40.0 


1 10" 


-10.3 


-9.3 


-3.9 


5.03 


1.50 


0.30 


50.0 


1103 


-8.2 


-8.1 


-4.0 


15.46 


19.72 


1.28 


50.0 


5103 


-8.8 


-8.4 


-4.0 


11.22 


13.68 


1.22 


50.0 


1 10" 


-9.1 


-8.4 


-4.0 


10.99 


11.84 


1.08 


50.0 


5 10" 


-9.6 


-8.5 


-4.0 


12.46 


8.00 


0.64 


50.0 


1105 


-9.7 


-8.7 


-4.0 


12.57 


6.30 


0.50 


50.0 


5 105 


-10.1 


-9.1 


-3.9 


8.65 


2.92 


0.34 


50.0 


1 10" 


-10.2 


-9.2 


-3.9 


6.33 


1.95 


0.31 


Go = 100 


2.0 


1 103 


-11.2 


-13.6 


-9.3 


1.08 


536.12 


498.10 


2.0 


5103 


-11.4 


-12.4 


-8.4 


2.97 


69.85 


23.49 


2.0 


1 10" 


-11.4 


-11.9 


-7.8 


3.43 


23.10 


6.74 


2.0 


5 10" 


-11.5 


-10.9 


-6.5 


2.10 


1.76 


0.84 


2.0 


1105 


-11.0 


-10.6 


-6.0 


1.21 


1.66 


1.37 


2.0 


5105 


-10.7 


-10.0 


-4.9 


0.28 


0.16 


0.56 


2.0 


1 10" 


-11.0 


-9.8 


-4.6 


0.22 


0.03 


0.13 


4.0 


1103 


-11.1 


-12.1 


-7.2 


0.29 


6.97 


23.75 


4.0 


5103 


-11.3 


-11.1 


-5.9 


0.18 


0.24 


1.31 


4.0 


1 10" 


-11.4 


-10.7 


-5.3 


0.13 


0.07 


0.58 


4.0 


510" 


-11.3 


-10.2 


-4.6 


0.14 


0.03 


0.20 


4.0 


1105 


-11.1 


-10.0 


-4.4 


0.18 


0.04 


0.23 


4.0 


5105 


-10.9 


-9.6 


-4.2 


0.28 


0.03 


0.11 
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Table E.l. Continued 







A'(HCO+)° 


X(HCN)° 


XCCO)" 


HCN/CO 


HCO+/CO 


HCO+/HCN 


(mag.) 


(cm-3) 








% 


% 




4.0 


1 10" 


-11.2 


-9.5 


-4.2 


0.32 


0.01 


0.04 


6.0 


110' 


-11.0 


-11.2 


-6.0 


0.15 


0.57 


3.73 


6.0 


510^ 


-11.1 


-10.6 


-4.9 


0.10 


0.07 


0.68 


6.0 


1 10" 


-11.1 


-10.3 


-4.6 


0.13 


0.05 


0.43 


6.0 


510* 


-11.2 


-10.0 


-4.3 


0.24 


0.04 


0.17 


6.0 


1105 


-11.1 


-9.8 


-4.2 


0.29 


0.05 


0.17 


6.0 


5105 


-11.0 


-9.6 


-4.1 


0.41 


0.04 


0.09 


6.0 


1 10* 


-11.3 


-9.5 


-4.1 


0.44 


0.02 


0.04 


10.0 


110^ 


-10.3 


-10.1 


-5.3 


0.92 


1.27 


1.39 


10.0 


510' 


-10.0 


-9.8 


-4.4 


0.72 


0.93 


1.30 


10.0 


1 10* 


-10.1 


-9.8 


-4.3 


0.64 


0.74 


1.16 


10.0 


510* 


-10.4 


-9.7 


-4.1 


0.59 


0.37 


0.63 


10.0 


5 10^ 


-10.8 


-9.6 


-4.0 


0.58 


0.12 


0.21 


10.0 


1 10* 


-11.0 


-9.6 


-4.0 


0.56 


0.07 


0.12 


16.0 


110' 


-8.8 


-8.5 


-4.3 


13.03 


14.15 


1.08 


16.0 


510^ 


-9.2 


-8.9 


-4.2 


5.81 


5.81 


1.00 


16.0 


1 10* 


-9.4 


-8.9 


-4.1 


4.71 


3.81 


0.81 


16.0 


510* 


-9.9 


-9.1 


-4.0 


3.32 


1.64 


0.50 


16.0 


110^ 


-10.0 


-9.2 


-4.0 


2.62 


1.19 


0.45 


16.0 


5 10^ 


-10.4 


-9.5 


-4.0 


1.35 


0.49 


0.36 


16.0 


1 10* 


-10.5 


-9.5 


-4.0 


1.04 


0.30 


0.29 


26.0 


1103 


-8.4 


-8.2 


-4.1 


15.16 


17.18 


1.13 


26.0 


510' 


-9.0 


-8.5 


-4.1 


9.34 


9.75 


1.04 


26.0 


1 10* 


-9.2 


-8.5 


-4.0 


8.40 


7.31 


0.87 


26.0 


510* 


-9.7 


-8.7 


-4.0 


7.50 


3.84 


0.51 


26.0 


1 10' 


-9.8 


-8.9 


-4.0 


6.53 


2.87 


0.44 


26.0 


5105 


-10.2 


-9.2 


-4.0 


3.42 


1.21 


0.35 


26.0 


1 10* 


-10.3 


-9.4 


-4.0 


2.42 


0.77 


0.32 


40.0 


1103 


-8.3 


-8.1 


-4.1 


15.43 


18.82 


1.22 


40.0 


5103 


-8.9 


-8.4 


-4.0 


10.58 


12.39 


1.17 


40.0 


1 10* 


-9.1 


-8.4 


-4.0 


9.92 


10.21 


1.03 


40.0 


510* 


-9.6 


-8.6 


-4.0 


10.22 


6.38 


0.62 


40.0 


1105 


-9.7 


-8.8 


-4.0 


9.69 


4.92 


0.51 


40.0 


5105 


-10.1 


-9.1 


-4.0 


5.80 


2.14 


0.37 


40.0 


1 10* 


-10.2 


-9.3 


-3.9 


4.12 


1.37 


0.33 


50.0 


1103 


-8.3 


-8.1 


-4.0 


15.49 


19.71 


1.27 


50.0 


5103 


-8.8 


-8.4 


-4.0 


11.06 


13.65 


1.23 


50.0 


1 10* 


-9.1 


-8.4 


-4.0 


10.55 


11.67 


1.11 


50.0 


510* 


-9.5 


-8.6 


-4.0 


11.47 


7.93 


0.69 


50.0 


1105 


-9.7 


-8.7 


-4.0 


11.28 


6.23 


0.55 


50.0 


5105 


-10.1 


-9.1 


-3.9 


7.26 


2.76 


0.38 


50.0 


1 10* 


-10.2 


-9.3 


-3.9 


5.23 


1.78 


0.34 



Notes. logarithmic relative abundance X(mol)= log(yV(mol)/A'(H2)) 
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Table E.2. Results of the PDR models with solar metallicities. 







X(HCO+r 


XCHCN)" 


X(COf 


HCN/CO 


HC07C0 


HCO+/HCN 


(mag.) 


(cm-3) 








% 


% 




Go = 10 


9 n 

z,.u 


1 10^ 


-116 
- 1 i .u 


-1 1 n 

- 1 1 .u 


-fl s 

u — J 


fl 84 
u.ot 


fl 4S 


n S4 




'i 10' 

J lU 


-1 1 fi 
- 1 i .u 


-in 4 


-S 1 


fl 1 S 

U. 1 J 


fl fl7 


0.12 


2.0 


1 10* 


-1 1.5 


-10.2 


-4.6 


0.1 1 


0.01 


0.10 


2.0 


5 10* 


-1 1.3 


-9.6 


-3.8 


0.18 


0.01 


0.07 




I 10^ 


-11.1 


-0 3 

y .3 


-3 7 


fl 7Q 


fl fl7 


n nfi 


2.0 


'i 10^ 


-1 1.2 


-8 Q 
o.y 


-3 S 

D.J 


fl SO 


fl fli 
u.u i 


n n7 


7 


1 10^ 

i 1\J 


-1 1.2 


-8 8 
o.o 


-3 4 
j.^ 


fl fiS 
u.uo 


fl fli 
u.u i 


n 01 


4 n 

t.o 


1 in3 

1 lU 


1 1 4 


1 n s 

lU. J 


4 7 


n 17 

u. iz, 


fl fl^ 
u.u J 


n 70 




103 


-11.1 


y .y 


-3 8 

J.O 


fl 10 

U. 1 7 


fl fl3 
u.u J 


n 16 

U. lU 


4.0 


1 10* 


-1 1.0 


-9.7 


-3.7 


0.29 


0.03 


0.12 


4.0 


5 10* 


-11.1 


-9.3 


-3.5 


0.64 


0.03 


0.05 


4. 


I 10^ 


-11.1 


-0 1 

y . 1 


-3 S 

J.J 


fl 83 

U.O J 


fl fl3 
u.u J 


n n4 


4. 


'i 10^ 


-113 
- 1 i.j 


-8 Q 


-3 4 
j.^ 


1.22 


fl fl7 


n ni 


4 fl 


i lU 


1 1 

- L L.J 


8 8 
o.o 


^ 4 


1 

1 .J J 


fl fli 
u.u i 


n ni 


u.u 


I 103 

1 J-U 


-ins 


-in 1 

- lU. 1 


-4.4 


fl 3n 

U. JU 


fl IS 

U. 1 J 


n SI 

U.J A 




103 

J lU 


-ins 

lU. J 


-Q 7 


-3 7 


fl 48 


fl 16 

U. lU 


n 34 

U. Jt 


fi n 

u.u 


1 Ifl-* 

i lU 


-ins 

lU. J 


-Q S 
y .J 


-3 fi 
^.u 


fl fiS 
u.u J 


fl 1 S 

U. 1 J 


n 93 


6.0 


5 10* 


-10.7 


-9.2 


-3.5 


1.03 


0.10 


0.10 


6.0 


1 10^ 


-10.7 


-9.1 


-3.4 


1.19 


0.09 


0.08 


fl n 

u.u 


105 

J J-U 


-1 1 n 

- 1 i .u 


-Q fl 

y .\j 


-3 4 

j.^ 


1 4fi 
1 .tu 


fl fl4 


n n3 

U.UJ 


fl fl 
u.u 


i lU 


1 1 1 


Q fl 

7.U 


^ 4 


1 S7 


fl fl^ 
U.U J 


n vo 


iU.U 


1 1fl3 

1 iU 


-Q 4 


-Q 7 


-3 7 


2.81 


3 07 


1 40 


1fl fl 

iU.U 


5 1fl3 

J iU 


-Q S 
y.j 


-Q 3 


-3 S 

J.J 


1.74 


2.12 


1 77 


Ifl fl 

iU.U 


1 Ifl** 

1 iU 


-Q 7 


-Q 3 


-3 S 


1 Sfi 
1 . JU 


1.44 


n 07 


10.0 


5 10* 


-10.1 


-9.2 


-3.4 


1.79 


0.64 


0.36 


10.0 


1 10^ 


-10.2 


-9.1 


-3.4 


1.99 


0.48 


0.24 


Ifl fl 

iU.U 


S 1fl5 

J iU 


-ins 

lU. J 


-Q fl 

7.U 


-3 4 


2.17 


fl 73 

U.Z. J 


n in 

U. lU 


1fl fl 

iU.U 


1 iU 


-in 7 

- lU. / 


-Q fl 

y .\j 


-3 4 


7 10 


fl 16 

U. lU 


n n7 

u.u / 


1 fl fl 
i u.u 


1 10^ 

1 iU 


-8 fi 
o.u 


-8 4 


-3 S 

J.J 


Ifl S8 

iU. JO 


1 1 00 
1 i .yy 


1 13 

1 . 1 J 


1 fl fl 

iU.U 


5 10^ 

J iU 


-0 1 

y. i 


-8 fi 
o.u 


-3 4 


fi 70 


S 78 
J . / o 


n SS 

U.OJ 


16.0 


1 10* 


-9.3 


-8.5 


-3.4 


6.63 


4.18 


0.63 


16.0 


5 10* 


-9.7 


-8.7 


-3.4 


5.79 


2.05 


0.35 


1 fi fl 
i u.u 


1 10^ 

1 iU 


-Q Q 
y.y 


-8 8 
o.o 


-3 4 


S fli 

J .U i 


1 SS 

1 .J J 


n 31 

U.J A 


16.0 


5 10^ 

J iU 


-10.2 


-Q fl 

7.U 


-3 3 
J.J 


3 07 

J .y 1 


0.72 


n 18 

U. lO 


1 fl fl 
i u.u 


1 10^ 

1 iU 


-in 4 


-8 7 
o. / 


-3 3 
J.J 


4 fl3 

T-.U J 


fl 40 


n 17 

U. Ix. 


Ifl fl 
zu.u 


1 1fl3 
1 iU 


8 4 


8 7 


^ 4 

J.H 


1 7 4fi 


14 40 


1 16 


26.0 


5 10^ 


-8.9 


-8.3 


-3.4 


9.03 


8.66 


0.96 


26.0 


1 10* 


-9.2 


-8.3 


-3.4 


9.03 


7.03 


0.78 


26.0 


5 10* 


-9.6 


-8.5 


-3.4 


8.94 


4.13 


0.46 


7fi fl 
.1.U.U 


1 10^ 

1 iU 


-Q 7 

y. 1 


-8 fi 
o.u 


-3 4 


8 76 
o.zu 


3 71 


n 30 


7fi fl 


5 10^ 

J iU 


-in 1 

- lU. i 


-8 Q 
o.y 


-3 3 
J.J 


fi 74 

U. /M' 


1 S3 
1 .J J 


n 93 


7fi fl 


1 10^ 

1 iU 


-10.2 


-8 7 
o. / 


-3 3 
J.J 


fi 86 
u.ou 


1 flS 

1 .UJ 


n IS 

u.u 


4fl fl 


1 iU 


-8 3 


-8 1 

0. 1 


-3 4 
j.t 


1 3 1 S 
i J. 1 J 


16.22 


1 73 


4fl fl 


S 1fl3 
^ iU 


-8 Q 

Q.y 


-8 7 


-3 4 
j.t 


10.14 


in 78 

lU. / o 


1 n6 

1 .uu 


4fl fl 


1 10^ 

1 iU 


-Q 1 


-8 7 


-3 4 
j.t 


Ifl 38 

iU. JO 


44 


n 01 


40.0 


5 10* 


-9.5 


-8.4 


-3.3 


11.17 


6.46 


0.58 


40.0 


1105 


-9.7 


-8.6 


-3.3 


10.92 


5.20 


0.48 


40.0 


5105 


-10.0 


-8.9 


-3.3 


9.75 


2.58 


0.27 


40.0 


1 10" 


-10.2 


-8.7 


-3.3 


10.02 


1.78 


0.18 


50.0 


1103 


-8.3 


-8.1 


-3.4 


13.41 


17.13 


1.28 


50.0 


510' 


-8.8 


-8.2 


-3.4 


10.62 


11.84 


1.11 


50.0 


1 10* 


-9.0 


-8.2 


-3.4 


11.01 


10.68 


0.97 


50.0 


510* 


-9.5 


-8.4 


-3.3 


12.26 


7.85 


0.64 


50.0 


1105 


-9.6 


-8.6 


-3.3 


12.26 


6.44 


0.53 


50.0 


5105 


-10.0 


-8.8 


-3.3 


11.49 


3.28 


0.29 


50.0 


1 10" 


-10.1 


-8.7 


-3.3 


11.88 


2.28 


0.19 


Go = 50 


2.0 


110^ 


-11.4 


-11.9 


-8.0 


2.55 


17.66 


6.93 


2.0 


5103 


-11.7 


-11.0 


-6.6 


1.14 


0.47 


0.41 


2.0 


1 10* 


-11.7 


-10.7 


-6.0 


0.68 


0.15 


0.23 


2.0 


510* 


-11.5 


-10.1 


-4.7 


0.18 


0.03 


0.14 


2.0 


1105 


-11.4 


-9.9 


-4.3 


0.12 


0.01 


0.12 


2.0 


5105 


-10.8 


-9.2 


-3.7 


0.23 


0.02 


0.07 
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Table E.2. Continued 







X(HCO+)'' 


X(HCN)° 


XCCO)" 


HCN/CO 


HC07C0 


HCO+/HCN 


(mag.) 


(cm-3) 








% 


% 




2.0 


1 10" 


-9.4 


-7.1 


-3.6 


0.85 


0.02 


0.02 


4.0 


1103 


-11.5 


-11.0 


-5.7 


0.14 


0.08 


0.58 


4.0 


5103 


-11.5 


-10.3 


-4.3 


0.08 


0.01 


0.14 


4.0 


1 10* 


-11.4 


-10.1 


-4.0 


0.12 


0.01 


0.11 


4.0 


510* 


-11.2 


-9.5 


-3.7 


0.32 


0.02 


0.06 


4.0 


1105 


-11.2 


-9.4 


-3.6 


0.43 


0.02 


0.04 


4.0 


5105 


-11.0 


-9.0 


-3.5 


0.73 


0.02 


0.02 


4.0 


1 10" 


-10.2 


-7.5 


-3.5 


1.16 


0.01 


0.01 


6.0 


1103 


-11.2 


-10.5 


-4.9 


0.17 


0.07 


0.45 


6.0 


5103 


-10.9 


-10.0 


-3.8 


0.23 


0.06 


0.27 


6.0 


1 10* 


-10.9 


-9.8 


-3.7 


0.33 


0.06 


0.18 


6.0 


510* 


-11.0 


-9.4 


-3.5 


0.62 


0.05 


0.08 


6.0 


1105 


-11.0 


-9.3 


-3.5 


0.76 


0.04 


0.06 


6.0 


5105 


-11.0 


-9.0 


-3.4 


1.02 


0.03 


0.03 


6.0 


1 10" 


-10.2 


-7.9 


-3.4 


1.25 


0.02 


0.02 


10.0 


1103 


-10.1 


-9.7 


-4.0 


1.11 


0.85 


0.77 


10.0 


5103 


-9.9 


-9.5 


-3.6 


1.02 


0.95 


0.93 


10.0 


1 10* 


-9.9 


-9.5 


-3.5 


0.97 


0.73 


0.75 


10.0 


510* 


-10.3 


-9.3 


-3.4 


1.05 


0.36 


0.34 


10.0 


1105 


-10.4 


-9.3 


-3.4 


1.12 


0.27 


0.24 


10.0 


5105 


-10.7 


-9.2 


-3.4 


1.27 


0.12 


0.10 


10.0 


1 10" 


-10.5 


-8.0 


-3.4 


1.55 


0.08 


0.05 


16.0 


1103 


-8.8 


-8.6 


-3.6 


9.36 


10.90 


1.16 


16.0 


5103 


-9.2 


-8.8 


-3.5 


4.95 


4.54 


0.92 


16.0 


1 10* 


-9.4 


-8.7 


-3.4 


4.67 


3.21 


0.69 


16.0 


510* 


-9.8 


-8.9 


-3.4 


3.76 


1.58 


0.42 


16.0 


1105 


-9.9 


-9.0 


-3.4 


3.20 


1.19 


0.37 


16.0 


5105 


-10.2 


-9.0 


-3.4 


2.60 


0.53 


0.20 


16.0 


1 10" 


-10.2 


-8.0 


-3.4 


2.98 


0.34 


0.11 


26.0 


1103 


-8.5 


-8.3 


-3.5 


12.16 


14.10 


1.16 


26.0 


5103 


-9.0 


-8.4 


-3.4 


7.76 


7.55 


0.97 


26.0 


1 10* 


-9.2 


-8.4 


-3.4 


7.55 


6.07 


0.80 


26.0 


510* 


-9.6 


-8.6 


-3.4 


6.94 


3.61 


0.52 


26.0 


1105 


-9.7 


-8.7 


-3.4 


6.15 


2.81 


0.46 


26.0 


5105 


-10.1 


-8.9 


-3.4 


4.94 


1.26 


0.26 


26.0 


1 10" 


-10.2 


-8.3 


-3.3 


5.29 


0.81 


0.15 


40.0 


1103 


-8.3 


-8.2 


-3.4 


12.97 


15.95 


1.23 


40.0 


5103 


-8.9 


-8.3 


-3.4 


8.88 


9.64 


1.08 


40.0 


1 10* 


-9.1 


-8.3 


-3.4 


8.88 


8.42 


0.95 


40.0 


510* 


-9.5 


-8.5 


-3.4 


9.11 


5.91 


0.65 


40.0 


1105 


-9.6 


-8.7 


-3.4 


8.64 


4.76 


0.55 


40.0 


5105 


-10.0 


-8.9 


-3.3 


7.58 


2.20 


0.29 


40.0 


1 10" 


-10.1 


-8.3 


-3.3 


8.18 


1.43 


0.17 


50.0 


1103 


-8.3 


-8.1 


-3.4 


13.25 


16.90 


1.28 


50.0 


5103 


-8.8 


-8.3 


-3.4 


9.35 


10.66 


1.14 


50.0 


1 10* 


-9.0 


-8.3 


-3.4 


9.47 


9.61 


1.01 


50.0 


510* 


-9.5 


-8.5 


-3.4 


10.14 


7.29 


0.72 


50.0 


1105 


-9.6 


-8.6 


-3.3 


9.92 


5.99 


0.60 


50.0 


5 105 


-10.0 


-8.9 


-3.3 


9.16 


2.84 


0.31 


50.0 


1 10" 


-10.1 


-8.5 


-3.3 


9.77 


1.85 


0.19 


Go = 100 


2.0 


1103 


-11.5 


-12.6 


-8.6 


2.27 


69.47 


30.54 


2.0 


5103 


-11.7 


-11.5 


-7.4 


1.97 


2.80 


1.42 


2.0 


1 10* 


-11.8 


-11.1 


-6.8 


1.46 


0.85 


0.59 


2.0 


510* 


-11.7 


-10.5 


-5.5 


0.47 


0.11 


0.23 


2.0 


1105 


-11.6 


-10.3 


-5.0 


0.23 


0.04 


0.19 


2.0 


5 105 


-10.5 


-9.5 


-4.0 


0.12 


0.05 


0.38 


2.0 


1 10" 


-10.7 


-9.2 


-3.8 


0.15 


0.01 


0.07 


4.0 


1103 


-11.6 


-11.3 


-6.3 


0.24 


0.27 


1.13 


4.0 


5103 


-11.7 


-10.6 


-4.7 


0.07 


0.01 


0.18 


4.0 


1 10* 


-11.6 


-10.3 


-4.3 


0.08 


0.01 


0.12 


4.0 


510* 


-11.4 


-9.7 


-3.8 


0.21 


0.01 


0.07 


4.0 


1105 


-11.3 


-9.5 


-3.7 


0.30 


0.01 


0.05 


4.0 


5 105 


-10.7 


-9.1 


-3.6 


0.53 


0.03 


0.06 
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Table E.2. Continued 







X(HCO+)'' 


X(HCN)" 


XCCO)" 


HCN/CO 


HC07C0 


HCO+/HCN 


(mag.) 










% 


% 




4.0 


1 10^ 


-11.0 


-9.0 


-3.5 


0.60 


0.01 


0.02 


6.0 


1103 


-11.3 


-10.7 


-5.1 


0.13 


0.06 


0.47 


6.0 


510' 


-11.1 


-10.1 


-3.9 


0.16 


0.04 


0.26 


6.0 


1 10* 


-11.0 


-9.9 


-3.8 


0.23 


0.04 


0.18 


6.0 


510* 


-11.1 


-9.5 


-3.6 


0.45 


0.04 


0.08 


6.0 


110^ 


-11.1 


-9.4 


-3.5 


0.56 


0.03 


0.06 


6.0 


5105 


-10.8 


-9.1 


-3.5 


0.80 


0.04 


0.05 


6.0 


1 10** 


-11.1 


-9.0 


-3.4 


0.83 


0.02 


0.02 


10.0 


110^ 


-10.3 


-9.8 


-4.2 


0.90 


0.61 


0.68 


10.0 


510' 


-10.0 


-9.6 


-3.6 


0.83 


0.64 


0.76 


10.0 


1 10* 


-10.1 


-9.5 


-3.6 


0.82 


0.52 


0.64 


10.0 


510* 


-10.3 


-9.4 


-3.5 


0.88 


0.28 


0.32 


10.0 


110' 


-10.5 


-9.4 


-3.4 


0.93 


0.21 


0.23 


10.0 


510^ 


-10.7 


-9.2 


-3.4 


1.01 


0.11 


0.11 


10.0 


1 10** 


-10.8 


-9.2 


-3.4 


0.99 


0.06 


0.06 


16.0 


110^ 


-8.8 


-8.7 


-3.7 


8.80 


10.47 


1.19 


16.0 


5103 


-9.2 


-8.9 


-3.5 


4.15 


4.05 


0.98 


16.0 


1 10* 


-9.4 


-8.8 


-3.5 


3.85 


2.85 


0.74 


16.0 


510* 


-9.8 


-9.0 


-3.4 


3.03 


1.40 


0.46 


16.0 


110' 


-10.0 


-9.0 


-3.4 


2.59 


1.06 


0.41 


16.0 


510^ 


-10.3 


-9.1 


-3.4 


2.15 


0.47 


0.22 


16.0 


1 10** 


-10.4 


-9.1 


-3.4 


2.13 


0.29 


0.14 


26.0 


110^ 


-8.5 


-8.3 


-3.5 


12.27 


14.21 


1.16 


26.0 


5103 


-9.0 


-8.4 


-3.4 


7.31 


7.17 


0.98 


26.0 


1 10* 


-9.2 


-8.4 


-3.4 


7.04 


5.74 


0.82 


26.0 


510* 


-9.6 


-8.6 


-3.4 


6.23 


3.45 


0.55 


26.0 


1105 


-9.7 


-8.8 


-3.4 


5.44 


2.68 


0.49 


26.0 


5105 


-10.0 


-9.0 


-3.4 


4.41 


1.18 


0.27 


26.0 


1 10** 


-10.2 


-8.9 


-3.4 


4.47 


0.73 


0.16 


40.0 


1103 


-8.3 


-8.2 


-3.4 


13.13 


16.13 


1.23 


40.0 


5103 


-8.9 


-8.3 


-3.4 


8.45 


9.27 


1.10 


40.0 


1 10* 


-9.1 


-8.3 


-3.4 


8.37 


8.09 


0.97 


40.0 


510* 


-9.5 


-8.5 


-3.4 


8.38 


5.76 


0.69 


40.0 


1105 


-9.6 


-8.7 


-3.4 


7.86 


4.64 


0.59 


40.0 


5105 


-10.0 


-8.9 


-3.3 


6.96 


2.10 


0.30 


40.0 


1 10** 


-10.1 


-8.9 


-3.3 


7.22 


1.32 


0.18 


50.0 


1103 


-8.3 


-8.2 


-3.4 


13.42 


17.11 


1.27 


50.0 


5103 


-8.8 


-8.3 


-3.4 


8.91 


10.28 


1.15 


50.0 


1 10* 


-9.0 


-8.3 


-3.4 


8.94 


9.28 


1.04 


50.0 


510* 


-9.4 


-8.5 


-3.4 


9.40 


7.14 


0.76 


50.0 


1105 


-9.6 


-8.7 


-3.4 


9.10 


5.88 


0.65 


50.0 


5 105 


-9.9 


-8.9 


-3.3 


8.51 


2.71 


0.32 


50.0 


1 10^ 


-10.1 


-8.8 


-3.3 


8.92 


1.72 


0.19 



Notes. logarithmic relative abundance X(mol)= log(A'(mol)/N(H2)) 



